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1. Introduction. 

1.1. Motivation and goal. This article is concerned with the existence and 
uniqueness of global regular solutions to a classical 3-D thermoviscoelastic system 
at small strains. The system describes materials which have the properties both of 
elasticity and viscosity. Such materials are usually referred to as Kelvin- Voigt type. 

As noted in the recent paper on this subject by Roubi'cek [21] - and according to 
our best knowledge as well - the existence of global solutions to a thermoviscoelastic 
system with constant both specific heat and heat conductivity is, in spite of great 
effort through many decades, still open in dimensions n > 2. In dimension n = 1 it 
was established in the pioneering papers by Slemrod [22] , Dafermos [5] and Defermos- 
Hsiao [6]. 

The local in time existence and global uniqueness of a weak solution to 3-D ther- 
moviscoelastic system with constant specific heat and heat conductivity has been 
proved by Bonetti-Bonfanti [3]. Other known results on multidimensional thermo- 
viscoelasticity deal with a modified energy equation. Modifications involve either 
nonconstant specific heat or nonconstant heat conductivity. A thermoviscoelastic 
system with temperature-dependent specific heat has been addressed by Blanchard- 
Guibe [2] where the existence of global, weak-renormalized solutions has been proved, 
and recently in [21] where the existence of a very weak solution has been established. 
We mention also that the framework of renormalized solutions has been applied in 
[26] for 3-D thermoviscoelastic system arising in structural phase transitions. 

In a more general setting allowing for large strains a 3-D thermoviscoelastic sys- 
tem has been studied under small data assumption by Shibata [23] and recently by 
Gawinccki-Zajaczkowski [11]. 

For thermoviscoelastic problems with a modified heat conductivity we refer to 
Eck-Jarusek-Krbec [8] and the references therein. 
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In the present paper we consider a thermoviscoelastic system with specific heat 
linearly increasing with temperature and with constant heat conductivity. Such set- 
ting is a particular case of systems addressed in [2] and [21]. 

The novelty of the existence result presented in this paper concerns the regularity 
of a 3-D global solution corresponding to sufficiently smooth but arbitrary in size 
initial data. The proof of the existence theorem is based on the successive approxima- 
tion method. The key regularity estimates are derived with the help of the parabolic 
theory in anisotropic Sobolev spaces W^ {Sl T ), Q T =!lx (0,T), p,p € (l,oo), 
with a mixed norm with respect to space and time variables. Such framework has 
been previously applied by the authors [19] to the thermoviscoelastic system arising 
in shape memory alloys. It allowed to generalize the former results on this subject in 
[27]. ' 

As known, in deriving a priori estimates for a solution of a system of balance 
laws it is common to begin with estimates arising from the conservation of a total 
energy. Such estimates provide Loo-time regularity for the conserved quantities. To 
take advantage of such time regularity in deriving subsequent regularity estimates it is 
desirable to work in Sobolev spaces with a mixed norm, for example Wp 'p (Sl T ), where 
the space exponent p is determined by the energy structure and the time exponent po 
may be arbitrarily large. This is the idea behind using the framework of Sobolev spaces 
with a mixed norm to the thermoviscoelastic system under considerations. The theory 
of IBVP's in Sobolev spaces with a mixed norm is the subject of recent theoretical 
studies. We apply the general results due to Krylov [13] and Denk-Hieber-Priiss [7]. 

1.2. Thermoviscoelastic system. The system under consideration has the fol- 
lowing form: 

(1.1) u tt - V- [A lSt + A 2 (e - Oa)] = b, 

(1.2) c v d9 t - kA6 = -9(A 2 a) ■ e t + (Axet) ■ e t + g in n T = SI x (0, T), 
where 

e = e(«) = J(V«+ (Vm) t ), e t = e{u t ) = ^{Vu t + (V« t ) T ). 

Here O C M 3 is a bounded domain occupied by a body in a fixed reference con- 
figuration, and (0,T) is the time interval. The system is completed by appropriate 
boundary and initial conditions. Here we assume 

(1.3) u = 0, n-V6> = on S T = S x (0,T), 

(1.4) u\ t= o = u , u t \ t=0 =u 1 , 9\ t=Q = 9 Q in 0, 

where S is the boundary of ft and n is the unit outward normal to S. 

The field u : S1 T — > R 3 is the displacement, 9 : S1 T — > R + = (0, oo) is the absolute 
temperature, the second order tensors e — (eij)i,j=i,2,3 and et = {(£t)ij)i,j=i,2,3 
denote respectively the linearized strain and the strain rate. 

Equation (1.1) is the linear momentum balance with the stress tensor given by 
a linear thermoviscoelastic law of the Kelvin- Voigt type (cf. [8], Chap. 5.4) 



S = A l e t + A 2 {e-9cx). 
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The fourth order tensors A 1 = ((Ai)ijki)ij t k,l=i, 2,3 and 

A-2 = {(A2)ijki)i,j,k,l=i,2,3 are respectively the linear viscosity and the elasticity ten- 
sors, defined by 

(1.5) £ h-> A p £ = Xptrel + 2fx p e, p = 1,2, 

where Ai, \i\ are the viscosity constants, and A2, [ii are the Lame constants, both Ai, 
fi\ and A2, [ii with the values within the elasticity range 

(1.6) V P >0, 3Ap + 2^ p >0, p = l,2; 

= (<5ij)i.j=i.2,3 is the identity tensor. 

The second order symmetric tensor a — (aij)i,j=i,2 with constant , represents 
the thermal expansion. The vector field b : 51 T — > M 3 is the external body force. 

Above and hereafer the summation convention over the repeated indices is used, 
vectors and tensors are denoted by bold letters, and the dot denotes the inner product 
of tensors, e.g. 

(Ae) ■ e = Aijkiekieij. 

Moreover, 

Ae = {Ai jk ie k i)i,j=i,2,3 and V • (Ae) = ( -^-{A ijk ie k i) ) 

\ C ' x j / i=l,2,3 

Equation (1.2) is the energy balance in which the linear Fourier law for the heat 
flux, q = —k\76 with constant heat conductivity k > 0, and temperature-dependent 
specific heat, c v 6 with c v > 0, have been adopted. The first two terms on the right- 
hand side of (1.2) represent heat sources created by the deformation of the material 
and by the viscosity. The field g : il T — > K is the external heat source. 

The boundary conditions in (1.3) mean that the body is fixed at the boundary S 
and thermally isolated. The initial conditions (1.4) prescribe displacement, velocity 
and temperature at t = 0. 

The system (1.1)— (1.2) can be derived by various arguments of thermodynamics, 
see e.g. [9, 3]. In Section 2 we summarize its thermodynamic basis. 

1.3. Linear elasticity and viscosity operators. Assumptions. For further 
analysis we formulate problem (1.1)-(1.4) in terms of the linear viscosity and elasticity 
operators, Q l and Q 2 , defined by 

(1.7) u h- > Q p u = V • (A p e(u)) = /i p Au + (A p + /i p )V(V • it), p=l,2, 

with domains D(Q p ) = H 2 (Vl) n Hl(Q). 
Then system (1.1), (1.2) takes the form 



u tt - QiU t = Q 2 u - V • (8A 2 a.) + b, 

c v 68 t - kA6 = -0(A 2 a) ■ e t + (A^t) -£ t +9 in T , 



with boundary and initial conditions (1.3), (1.4). 
Throughout we shall assume that 

(AI) £1 C R 3 is a bounded domain with the boundary S of class at least C 2 ; T > 
is an arbitrary finite number; 
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(A2) a = (aij)ij = i,2,3 is a second order symmetric tensor with constant a^; 

(A3) The fourth order tensors A\ and A 2 are defined by (1.5) with the coefficients 
V P = !' 2 ' satisfying (1.6). 

We list the implications of assumption (A3) which are used in further analysis. 
The condition (1.6) ensures the symmetry of tensors A p : 

(1-9) (A p )ijki = (A p )jiki = (A p )kuj, p = 1,2, 

and their coercivity and boundedncss 

(1-10) (v|e| 2 < {Ape) -e < a* p \e\ 2 , p=l,2, 

where 

a pit = min{3A p + 2pb p , 2^ p }, a* = max{3A p + 2^ p , 2/j p }. 

Moreover, (1.6) ensures the following properties of operators Q pl p = 1, 2: 

- Q p are strongly elliptic (property holding true under weaker assumption fj, p > 
0, A p + 2fi p > 0, (see [20], Sect. 7)) and satisfy the estimate [17], Lemma 3.2: 

(1.11) c p \\u\\ H 2 {n) <\\Q p u\\ L2{n) for ueD(Q p ), p = 1,2, 

with positive constants c p depending on fi. Since clearly, 

\\QpU\\ L2 (n) < c p \\u\\ H 2^ n) , 

it follows that the norms ||QpM|U 2 (n) an d H M Hir 2 (fi) are equivalent on D(Q p ); 

- the operators Q p are self-adjoint on D(Q p ): 

(1 12) ^ pU ' ^ i2(n) = ~^ Vtt ' Vw )m«) - ( A p + M P )(V ■ u, V • u) i2 (n) 
= (u, Q p v) L2{Q) for u,v £ D(Q p ); 

- the operators — Q p are positive on D(Q p ): 

{-QpU, u) = M P l|Vu||| 2(n) + (A p + ^p)|| V • u||i a(n) > 
for u £ r>(Q p ). 

Hence, there exist fractional powers Q^j 2 with the domains 
D{Q]I 2 ) = ffj(fi), satisfying 

^ (Qp 7 V Qp /2 ^)i 2 (J2) = (-Qp", v)i 2 (Q) = («, -Q p f )i 2 (n) 

for u,f G D(Q p ). 

Let us also notice that by (1.10) and the Korn inequality 

(1.15) d 1/2 ||«||Hi ( n) < l|e(u)|U 2( ii) for u6ffj(n), d > 0, 

it follows that 

(1 ll<3p /2 «lll 2 (o) - MpI|Vm||| 2(0) + (Ap + Mp )||V • u||i a(n) 

= (A p e(ti),e(u)) i2( n) > <vll e ( u )lll 2 (n) > <vd||u|lffi(n)- 

Thus, the norms ||Qy 2 M||i 2 (n) an d IMIjj 1 ^) are equivalent on D(Q p ^ 2 ). 
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1.4. Main result. 

Theorem A. (Existence) Let the assumptions (Al)-(AS) hold, S G C 2 , T > 
finite and 

u eW 2 12 (a), UiSBgfO), foe<(fi), 
g t -^00,12 

where 6_ is a constant. Then there exists a solution to problem (1.1)-(1.4) such that 
u e C{[0,T];W 2 12 {Q)), u t € W 2 { 2 {tt T ), 9 e W 6 2,1 (Q T ) and 0(t) > #exp(-c T) = 
0* > 0, where the positive constant Co depends on ai*, a^, |a|, c„. 
Moreover, the following estimates are satisfied 

\\ u \\c([0,T];W 2 12 (n)) < c ll u t|lw^i (aT ), 

IWIw^ 1 ^) + H 6 'IU 6 2 ' 1 (n r ) ^ ^( T > \\ u o\\wl 2 (n) + Ihillu^o.) 
+ IMI 

where ip is an increasing positive function of its arguments. 

Theorem B. (Uniqueness) Let us assume that tensors A p , p = 1,2, satisfy 
(1.10). Then any solution (u,6) to problem (l.l)-(H) satisfying 

e t eL 2 (0,T;L 3 (n)), 
(1.17) e€L 2 {Q,T;L 00 {9)), 6 t e i 2 (0,T; L 3 (fi)), 

< 61* < 0, 

is uniquely defined. 

Corollary 1.1. TTie regular solution in Theorem A is uniquely defined. 

1.5. Relation to other results. We comment on the connections of our result 
to the two other global existence results in three space dimensions. Firstly, we mention 
the result by Roubfcek [21] who proved the existence of a very weak solution to the 
thermoviscoelasticity system (1.1)— (1.2) involving monotone viscosity of a p-Laplacian 
type, (Ai£ t ) ■ £t ~ \st\ p , and the specific heat having (u> — l)-polynomial growth, 
c*(0) <~ c t ,0'" _1 . This result, based on the Galerkin method, was obtained for L 1 -data 
under the conditions p > 2, w > 1 and p > 1 + (in 3-D). In the case of linear 
viscosity, p = 2, the latter condition implies that u) > 3/2, that is the growth of the 
specific heat is greater than 1/2. 

Our result concerns the case p — 2 and to = 2. We have to restrict ourselves to 
the linear viscosity, p = 2, because the proof relies on the results by Krylov [13] and 
Solonnikov [24] on the solvability of the linear problem 

Utt - V • Aie(« t ) = / 

with the boundary and initial conditions (1.3), (1.4) (see Lemma 3.4). 
Concerning the specific heat growth exponent, u) — 1, it seems that after some addi- 
tional technical effort it would be possible to admit u < 2. However, in the case of 
a constant specific heat, i.e., oj = 1, we have been faced with a serious mathematical 
obstacle. 
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As already mentioned in Subsection 1.1, the local existence result in such a case was 
obtained by Bonetti and Bonfanti [3] . 

Secondly, we recall the multidimensional result by Blanchard and Guibe [2] who 
addressed problem (1.1)— (1.2) equally in the prototype case p = 2 and to = 2, and in 
the more general setting involving linear viscosity, p — 2, specific heat with {u> — 1)- 
polynomial growth and a nonlinear thermoelastic coupling; more precisely, the term 
V6* in (1.1) was replaced by V/(6*) with / hawing an a-polynomial growth. The 
existence of solutions in the weak-renormalized sense was proved there by the Schaudcr 
fixed point theorem. It is worth to remark that in the case of the linear thermoelastic 
coupling, a = 1, the result in [2] requires the specific heat to have growth of the order 
greater than 1/2, as in the result by Roubicek [21]. 

1.6. Outline. In Section 2 we present a thermodynamic basis of system (1.1), 
(1.2). Section 3 recalls basic results on the Sobolev spaces with a mixed norm and 
on the solvability of boundary-value problems for linear parabolic equations in such 
spaces. In Section 4 we derive a priori estimates for problem (1.1)— (1.4). The pro- 
cedure consists in a recursive improvement of the basic energy estimates. The main 
tool in this procedure are the results on the solvability of linear parabolic problems 
in Sobolev spaces with a mixed norm. Section 5 presents the proof of Theorem A, 
which is based on the successive approximation method. The proof of the uniqueness, 
stated in Theorem B, is given in Section 6. 

Since a priori estimates in Section 4 are crucial for the proof of the global existence 
we advertise here the main steps of the procedure of deriving such estimates. First 
we prove the energy type estimate (see Lemma 4.2) 

(1.18) IKIU 2lOC (fi-) + IHk^Q-) + 11^11^,00(0-) < data. 

In Lemma 4.6 we show the estimates 

ll M *IU 2 ,oo(o T ) + IMIz, co (o,T;.H- 1 (n)) + \\ u t\\L 2 (o,T;H 1 (tn)) < c ll#ll l 2 (o t ) + data, 

and 

IKIUooCo.Tjff 1 ^)) + IMIlooCctw 2 ^)) + \\ u t\\L 2 (o,T;H 2 (n)) < c l|V6>|| £2(n T) + data. 

The norms of 9 will be later removed by some interpolation inequalities based on 
estimate (1.18). 

In Lemma 4.7 we obtain the estimate 

(!-!9) \\0\\ Loo (o,T;L 3 (ci)) + ||0||L 2 (o,T;ffi(n)) + IMIv 2 (n r ) < dat a, 

and next in Lemma 4.8, 

l|0t||L 2( n-) + ||V^|| Loo(0 ,T;i 2 (n)) < data. 
To deduce the boundedncss of 9 we first prove in Lemma 4.10 the estimate 
(1-20) l|0|U r .oo(n*') < data, r < oo, 

and in Lemma 4.17, 



(1.21) 



l|0|Uo.(ii T ) < data. 
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To get (1.21) we make use of the important inequality (see Lemma 4.9) 
\\et\\i. p , a (nT) < c{t){W\\L p ,An T ) +data), p,a G (l,oo), 

and the fact that the coefficient near t in (1.2) is proportional to 9. 

To establish the continuity of 9 (proved in Lemma 4.10) we need (1.21), estimates 

Hetllz-^o.TjiooCn)) < data (see Corollary 4.16), and 

ll^ll w 2 2,1 (J2 T ) — data ( see Corollary 4.18). 

Having the previous estimates for e t and the continuity of we finally prove in Lemma 
4.23 that 

9 G W 6 2,1 (0 T ), u t G W^(fl r ). 

2. Thermodynamic basis. System (1.1), (1.2) represents balance laws for the 
linear momentum and energy in a referential description, with the referential mass 
density assumed constant, normalized to unity, p = 1: 

(2.1) u tt -V.S = b, 

e t + V • q- S ■ e t = g, 

where S is the stress tensor, q - the referential heat flux, and e - the specific internal 
energy. 

The system is governed by two thermodynamic potentials: the free energy / = 
f(e,9), which by a thermodynamic requirement is strictly concave with respect to 
6, and the dissipation potential (called pseudopotential of dissipation in [10], [3]) 
T> = T>{et, V0; s, 9), which by a thermodynamic requirement is nonncgative, convex 
in (e t , V0) and such that X>(0, 0; e, 9) = 0. 

In the case of (1.1), (1.2) the free energy is specified by 

(2.2) f(e,9) = .U(9) + W(e,9), 
where 

(2.3) /*(0) = -^c v 9 2 , Cy = const > 0, 
is the caloric energy, and 

W(e, 9) = he- Oct) ■ A 2 (e - 9a) - ^a ■ (A 2 a) 

(2.4) 2 2 

= -s • (A 2 e) - 9e ■ (A 2 a) 

is the elastic energy; we recall that A 2 stands for the fourth order elasticity tensor 
and a for the second order thermal expansion tensor. 

The caloric energy (2.3) is associated with temperature-dependent caloric specific 
heat 

(2.5) c49) = -9f':(0)=c v 9, 
which gives rise to the term c v 99 t in energy equation (1.2). 
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We remark that in the case /* is given by the standard formula 

(2.6) U{6) = -c v 6\og?- + c v e + c, 

where c v , 9\, c are positive constants, the caloric specific heat is constant 

(2.7) c* = -6f:{6) = c v . 

This gives rise to the usual parabolic term c v 9 t in (1.2) in place of c v 99t- As mentioned 
in Section 1, in the case of a constant caloric heat there are serious mathematical 
obstacles in the proof of the global existence. 

A thermoviscoelastic system with the specific heat c*(#) given by (2.5) has been 
considered in [9], [2] and [21], where also more general forms of c*(#) have been 
analysed. Moreover, we mention that a fourth order thermoviscoelastic systems with 
temperature-dependent specific heat, arising in shape memory materials, have been 
studied in [27] and [19]. 

The dissipation potential corresponding to system (1.1), (1.2) is given by 



1 , . , k 



(2.8) V= -e t -(A 1 e 1 ) + 



26 v ' 2 



2 



l £t .(A l£t ) + ^|Vlogtf| 2 , 



where A\ is the viscosity tensor and k > the constant heat conductivity. 

In accord with the basic thermodynamic relations the internal energy e and the 
entropy r\ are related to the free energy / by the equations 

(2-9) e = f + 6r), r, = -/,,. 

For the free energy / defined by (2.2)-(2.4) this gives 

(2.10) e = ic,„6» 2 + -e- (A 2 e), rj = c v 6 + (A 2 a) ■ e. 

As a consequence of the second law of thermodynamics expressed by the Clausius- 
Duhem inequality, the stress tensor S and the heat flux q satisfy the following rela- 
tions: 

(2.11, + 



de de t ' H dVy 

For / given by (2.2)-(2.4) and V by (2.8) the formulas (2.11) yield the standard 
forms of the stress tensor and the heat flux 

(2.12) S = A 2 {e - 0a) + A lSt , q = k9 2 v\ = -fcV '6 : 

9 

Thus, S consists of two terms: the nondissipative equilibrium term determined by 
/, and the dissipative one determined by V. The dissipative heat flux q is entirely 
determined by V. 

Inserting the relations (2.10)i and (2.12) into balance laws (2.1) one arrives at 
the system (1.1)-(1.2). 

For further purposes (see Lemma 4.2) it is of interest to notice that on account 
of the identity 

df 

e t = (f + 9 v ) t = f t + 9 tV + 9 Vt =9 Vt + -^-- e t , 
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along with the relation (2.11)i, the energy balance (2.1)2 admits the form 

dV 

(2.13) 6 m + V ■q = e—-£ t +g. 

OEt 

For D given by (2.8) this leads to the following equivalent form of equation (1.2): 

(2.14) 6r h -kA6={A 1 £ t )-£ t +g. 

Let us also notice that assuming 6 > and using (2. 1 1) 2 , the equation (2.13) may 
be expressed as 

(2-15) 11t + V -l =a+ % 



where 



dV _ 1 dV 



+ ^(A 1 e t )-e t >0 



is the specific entropy production. From (2.15) it follows that system (1.1), (1.2) 
complies with the Clausius-Duhem inequality 

(2-16) % + V-|>|. 

3. Notation and auxiliary results. 

3.1. Notation. Let SI C K" be an open bounded subset of R™, n > 1, with 
a smooth boundary S, and S] T =!!x (0,T), S T = S x (0,T), T > finite. 
We introduce the following spaces: W£ (SI), k e N U {0}, p G [l,oo) - the Sobolev 
space on SI with the finite norm 



1 /p 

I \ ^ / I nQ„.lp 1 



E / 



l«l<fef! 

where a — (a\, . . . , a„) is a multiindex, a* € N U {0}, |a| = a\ + a 2 + . . . + a n , 
D° = ^...a£; # fc (^) = ^ 2 fe (tt); L p , Po (0 T ) = L po (0,T;L p (tt)), p, Po € [l,oo) - 
the space of functions u : (0, T) — > L p (f2) with the finite norm 



l«llw« T ) = (/ ll«(t)ll? p(n) *) 



1/P0 



1/ 2 (^ T ) = J L oc (0,T;L 2 (fi)) n L 2 (0,T;H 1 (SI)) - the space of functions u : (0,T) 
H 1 (Sl) with the finite norm 

IMk(fi T ) = ess SU P ll"(*)IU 2 (n) + ||Vu|| i2(n T ); 
te[o,T] 

V 2 '°(Sl T ) = V 2 (Sl T ) n C([0,T];L 2 (SI)) - the space with the finite norm 
IMIv^n*') = 4 m a ^ ll«(*)IU a (n) + I|Vu|| L2(0 t ); 
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Wp,'p /2 (n T ), fc,fc/2eNU{0}, p,p e [1, oo) - the Sobolev space with a mixed norm, 
which is a completion of C°°(f2 T )-functions under the finite norm 



l 

J ( E / \ D td>\ p dx 



|a|+2a<fc n 



Po/p \ 1/po 

dt 



p,p e [l,oo) - the Sobolev-Slobodecki space with the finite 

norm 



l u llw p ";; /2 (or) 



+ 



= £ P*^lk, P0 (n-) 

|a|+2o<[s] 



E 

"f! f! |a|+2o=[s] 



T T 



+ 



oo n 



E 

|a|+2o=[s] 



|j?gdMM)-^M^)l p 

| x _ x '\n+p(s-[s]) 

\D2d?u(x,t) - D2d$u(x,t')\r 

lt-t'| 1+ Kf-[f]) 



dxdx' 



Po/p 



dt 



dx 



po/p 



dtdt' 



i/p 



1/po 



where a e NU {0} and [s] is the integer part of s. For s odd the last term in the above 

norm vanishes whereas for s even the two last terms vanish. 

B l (0), i G K+, p,po € [l,oo) the Bcsov space with the finite norm 



B l (12) 



Up(o) + ( E / 



h l+(l-k)p 



where: 

fceNU{0}, meN, m>Z-fc>0, 

A^(ft, fi)tt, j £ N, h € K+, is the finite difference of the order 
j of the function u(x) with respect to Xj, with 
A}(h,Q)u = Ai(h,n)u 

= u(xi, . . . , Xi-i,Xi + h, x i+ i, ...,x n )- u(xi, . . . , x n ), 

Ai(h,n)u = A i (ft,fi)Aj _1 (/i,fi)u, 

and 

A\ (h, Sl)u = for x + jh £ ft. 

In [12] it has been proved that the norms of the Besov space B ppo (fl) are equiv- 
alent for different m and k satisfying the condition m > I — k > 0. 

By c we denote a generic positive constant which changes its value from formula to 
formula and depends at most on the imbedding constants, constants of the considered 
problem and the regularity of the boundary. 

By ip = <p((Ti, . . . , Ufe), k G N, we denote a generic function which is a positive 
increasing function of its arguments o\ , . . . , <7fc, and may change its form from formula 
to formula. 
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3.2. Auxiliary results. We need the following interpolation lemma 
Lemma 3.1. [1, Chap. 4, Sect. 18] Let u <= Wp^ 2 ^), s G R +; p,p Q G [l,co], 
0, C M 3 . Let a G K+ U {0}, and 

3 2 3 2 
k= — I \- \a\ + 2a + a < s. 

p pa q qo 

Then D"d?u G Wq'qJ 2 (Q T ), q > p, qo > Po, and there exists e G (0, 1) such that 
\\ D x d >\\wZrt 2 (nT) < £S ~ K W u \\w;-X 2 (nT) + ce~ K \HL p , P0 (nn- 

We recall from [4] the trace and the inverse trace theorems for Sobolcv spaces 
with a mixed norm. 

Lemma 3.2. (Traces in Wp% l2 {VL T )) 

(i) Let u e w^ 2 ^), s e k +; p , Po e (i,oo). 

Then u(x,to) = u(x,t)\t=t f or *o € [0, T] ; belongs to Bp^ p 2 J Pa (O), and 



where constant c does not depend on u. 

(ii) For a groen u G Bp, p 2 J Po (Q), , 
a function u G Wp'pl 2 (fl T ) such that u\ t= t„ — u for t G [0,T], and 



(ii) For a gi-uen u G -Bp lP o (f2), s G s > 2/po, P,Po G (l>oo)> ^ere exists 



where constant c does not depend on u. 

We recall also (see [1]) that if I > 1/p then every function from B l ppo (VL) has a 

trace on the boundary S belonging to B l PlP 1 ' p (S), and 

We apply the following imbeddings between Besov spaces. 

Lemma 3.3. [25, Theorem 4-6-1] Let C W 1 be an arbitrary domain. 

(a) Let s G R+, s > 0, p G (l,oo) and 1 < q\ < q2 < co. TTien 

C B^fi) C B'^fi) C fl'^fi) C B' Pt00 (Q) C B^ £ (0). 

f&j Le£ oo > q > p > 1, 1 < r < oo, < t < s < oo, and 

n n 

t H <s. 

P q 



Then 



B' Pir (Q) C Bj ir (fi). 



We recall now from [19] a result on the solvability of a linear parabolic system 
with elasticity operator Q in Sobolev space with a mixed norm. This result will be 
repeatedly used in Section 4 in deriving a priori estimates for viscoelasticity system 
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(1.1). It generalizes the result by Krylov [13] from the single parabolic equation to 
the following parabolic system 

u t Qu = f in T = n x (0, T), 

(3.1) u = Q on S T = Sx (0,T), 

u\ t =o = u in fl, 

where SI C R 3 , S = dfl, f = (/,), and 

Qu = u.Au + tA7(V • u) 
with /i > 0, v > 0. Let us notice that letting 

Q = Q\ 1^ = 1^1, v = \i+^\, 

the assumption (1.6) implies that \i > and > 0. 

Lemma 3.4. (Parabolic system in W 2 p l (n T ) [13, 19, 24]) 

(i) Assume that f e i P ,p„(« T ), «o € B^/ po (0) ; p, Po e (l,oo), S E C 2 . 7/ 
2 — 2/po — 1/p > £/ie compatibility condition u \s = is assumed. Then there exists 
a unique solution to problem (3.1) such that u E W 2 'p (Q T ) and 

(3-2) \\u\\ w 2,i g(nT) < c(||/|| ippo(nT) + K|| B 2-2/ P0(n) ) 

with a constant c depending on fl, S, p, Pq. 

(ii) Assume that f = V ■ g + b, g = {gij), b = (b l ), g,b E L PtPa {VL T ), it E 
B 2 ~ p 2 J Po (£1). Assume the compatibility condition 

u \s = if 1 - 2/po ~ 1/P > 0. 
Then there exists a unique solution to (3.1) such that u E VFp^ 2 (0 T ) and 

IMI^w^t) < c{\\g\\ Lppo{ nT) + WHi^itiT) 

+ \\u \\ B 2-y PO{n) ) 

with a constant c depending on fl, S, p, Pq. 

In the proof of Theorem A we shall apply also the following regularity result for 
a linear parabolic equation. This result is the special case of a more general theorem 
due to Denk-Hieber-Priiss [7, Theorem 2.3]. 

Lemma 3.5. (Parabolic equation in W p ; pa (£l T )) Let us consider the problem 

6 t - qA6 = g in J1 T , 
(3.4) n-V6» = on S T , 

6\ t=0 = 9 Q in O, 

where g(x,t) is a continuous function on il T such that inf q g > 0. Assume that 
g E L p ^ Po (Q T ), 6 a E Bp^/ Pa (il), p,p E (l,oo), S E C 2 , and the corresponding 
compatibility condition is satisfied. Then there exists a unique solution to problem 
(3.4) such that 9 E W 2 ^ o (Q T ) and 

( 3 - 5 ) W \\w2;l o (nn < c(\\g\\ Lppo{nT) + ||0 o || s 2-2 /PO(a) ) 

with a constant c depending on Q, T , S , inf^r g and sup^r g. 

Remark 3.6. The constants c in Lemmas 3.4 and 3.5 do not depend on T . ForT 
small the proof of this fact is evident whereas for T large it can be deduced by applying 
the same arguments as in the proof of Theorem 3.1.1 in [26, Ch. 3]. 
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4. A priori estimates. In this section we derive a priori estimates for solutions 
of problem (1.1)— (1.4) on an arbitrary finite time interval (0,T). The estimates are 
essential for the existence proof by the successive approximation method, presented 
in Section 5. 

The procedure of deriving a priori estimates consists in a recursive improvement 
of the basic energy estimates. The main tool used in this procedure is Lemma 3.4 
which provides the solvability of the viscoelasticity system (1.8)i in the Sobolev space 
Wp'p (Q T ). The applied procedure is aimed to establish the continuity of temperature 
and finally to apply Lemma 3.5 on the solvability of parabolic equation in Sobolev 
space Wl> q \(Q T ). 

Throughout this section we assume that assumptions (A1)-(A3) (see Sect. 1.3) 
hold, and 

(4.1) O >6>O in SI, g>0 in Q T , 

where is a positive constant. 

First we prove the lower bound on by using similar arguments as in [16, Lemma 
3.7], [27, Lemma 3.3]. 

Lemma 4.1. (Lower bound on 0) Let us assume that (4-1) holds. Then there 
exists a positive constant c depending only on parameters ai*, (from (1.10), \a\ 
(see A2), c v (see 1.2)), such that 

(4.2) 6{t) > 0exp(-cT) = 0* > for t £ [0,T\. 



Proof. For m £ K+ let us define the truncation 

m = max • 

and 



fl m (t) = jz £ n : 0(x,t) > ^ j • 



Multiplying (1.2) by — m e with g > 2 (admissible test function) and integrating over 
tt m (t) gives 



(4.3) 



(Ai£ t ) • s t 

n m (t) n m (t) o m (t) 



e dx + k I 0~ e A0dx+ I v L e ; — -dx 



J ~e\ dx= J -^( A 2a) ■ £tdx. 



n m (t) o m (t) 
The first term on the left-hand side of (4.3) is equal to 



(4.4) 



n m (t) n m (t) 



- Cv f d t e 2 -edx - Cv d ( 2 ~ e dx 
~ g - 2 J ° l m g-2dt J m ' 
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because d t 6^ e = 0fora;€fi\ Q m {t) = {x £ n : 6 m {t) = £}. 
The second term on the left-hand side of (4.3) equals 



(4.5) 



k J 0-BAe m dx = kj 6 m s A6 m dx = ^ k3 i)2 J v(-^rr 



dx, 



because V9 m = V8 for x £ fi m (i) and V# TO = for x £ fi \ Q, m {t). On account of 
(1.10) the third term on the left-hand side of (4.3) is bounded from below by 



(4.6) 



and the fourth one by 



/ 



{Met) ■ £t 

dm 



dx > au 



J 



■dx, 



.(*) 



.(*) 



(4.7) 



> 0. 



n m (t) 



In view of the boundcdncss of A2 and a, the integral on the right-hand side of (4.3) 
is estimated as follows: 



(4.8) 



n m (t) 



o m (t) n m (t) 



Now, setting 5/2 = ai* and incorporating (4.4)-(4.8) into (4.3) we arrive at 



g-2~dt 



dx 



<c J 9 2 m e dx <cj 6 2 m e dx, 
n m (t) n 



where in the last inequality we used the fact that 6 m > in £1 Hence, by the Gronwall 
inequality, it follows that 



J 6 2 m e{t)dx< J O^^Ojdxexp 



c{g - 2) 



for t£ [0,T], 



n n 
that is, 

(4-9) llC(*)lk- 2 (n) < HC(0)IU e _ a (n) exp(cT) 

with a constant c independent of g and to. Letting g — > 00, (4.9) yields the bound 

Mt) >0m(O)cxp(-cT). 
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Further, letting m — > oo and noting that for sufficiently large to, 9 m (0) = 
= niax{#o, ^} > we conclude the bound (4.2). □ 

Lemma 4.2. (Energy estimates) Let us assume that (4-1) holds, 9 > 0, and 

uoeH 1 ^), Ui gl 2 {Q), e eL 2 (n), beL 2A (n T ), g ^L x {n T ). 

Then a sufficiently smooth solution (u,0) to (t.l)-(H) satisfies the estimate 



(4.10) 



|N'|U 2 ,oo(fi*) + ll £ lli2,oo(^) + ll^lUa.ooCfJ*) 

+ \\e- 1 ve\\ L2m + \\e-^ 2 e t/ \\ L2{nt) 

< c(||e(«o)IU 2 (n) + ll«ilU 2 (fi) + ll^olU 2 (n) 
+ ||b|U 2 , l( n*) + lblU l( o*) + 1) = cA , 

where t <T. 

Proof. Note that the positivity of 9 is ensured by Lemma 4.1. Multiplying (1.1) 
by u t , integrating over f2 and integrating by parts using boundary condition (1.3)i, 
gives 



J \u t \ 2 dx + j (Ai£ t ) ■ e t dx + J A 2 (e - 9a) ■ e t dx 



(4.11) 



= J b ■ u t dx. 
n 

Further, integrating (1.2) over fi and by parts using (1.3)2 yields 
(4.12) ° v d 



2 dt 



j 9 2 dx + J ( 



9(A 2 a) ■ e t dx — J (Ai£ t ) ■ e t dx = J gdx. 
n n n n 

Next, multiplying (1.2) by 1/9, integrating over and integrating by parts leads to 

- r (Ai£ t ) • £ t 



(4.13) j t J [c v 9 + (A 2 a) ■ e]dx -k j ^-dx - J ■ 
n an 



-dx = j ^ (/.c. 



Adding by sides (4.11) and (4.12) gives 



(4.14) 



dtj 



±c v 9 2 + ^(A 2 e) -£+\\u t \ 2 



dx 



= J(bu t + g)dx. 



Now, multiplying (4.13) by a positive constant /3, and subtracting by sides from (4.14), 
we get 

d_ 
~dt 



I 



l -c v 9 2 + \{A 2 e) • e + i| M4 | 2 - /3c v 9 - p(A 2 a) ■ e 



dx 



(4.15) 



k(3 



9 2 



■dx + f3 



(Ai£ t ) ■ e t 



dx 



I 



b ■ u t + I 1 - £ 



dx. 



Integrating (4.15) with respect to time, noting that by the boundedncss of A 2 and a, 
^c v 9 2 + \{A 2 e) ■ e + i|« t | 2 - (3c v 9 - (5{A 2 a) ■ e 

- \ c -"° 2 + J( A 2£) -£ + ^\u t \ 2 -c, 



1G 



I. Pawkow and W. M. Zajczkowski 



and using the coercivity of Ai (see (1.10)), we get 



l|0(*)lli 2( n) + \\m\hn) + ll«*(*)lli 2 (n) + We-'vewl 



\L 2 (Q) T" N"tWIIX 2 (f2) 

9 o\\l 2 (n) + 

+ l|Mt'lli 2 ,oc(n' )II^IIl 2 , 1 (o*) + Wgh^sv) + c 



2(0*) 



||^ 1/2 £ t '|l!, W )<c(ll^o|lL ff2) + ||£(«o)H! 2W 



l«ill! 2 (o)) 



for t <T, with a constant c depending only on parameters. 

Now, applying the Young inequality to the second term on the right-hand side of the 
above inequality yields (4.10). □ 

Remark 4.3. By integrating the identity (4- 14) with respect to time one can 
immediately conclude that 



(4.16) 



K'|U 2 ,oc(f2*) + ||e|U 2i00 (f2t) + ||0|U 2iOO (f2*) < cA , 



which is a part of inequality (4-10). 

On the contrary to (4-10) this inequality does not require the assumption 9 > 0. We 
point out that in deriving further estimates we shall use just the bounds in (4-16) loos- 
ing the information contained in the two dissipative terms of (4-10). This information 
may be of importance in the analysis of the long time behaviour of solutions. 

Remark 4.4. We complement Lemma 4-2 by some physical interpretations. In 
view o/(2.10)i, (2.12) and boundary conditions (1.3), identity (4-H) represents the 
balance equation for the total energy 



d_ 

dt 



/H 



Ut\ 



dx 



-(Sn) ut+n - q]dS = J (b ■ u t + g)dx. 

n 



On the other hand, in view of (2.10)2, (2.12) 2 and the boundary condition (1.3)2, 
identity (4-13) represents the balance equation for the entropy 



J, I 



n ■ —dS = / adx 



^-dx 



1 



with the entropy production 

v=^\V0\ 2 +Q(A l£t )-e t >O. 
Equation (4-15) represents the so-called availability identity 



d 
It 



\\u t \ 



flrj^J dx + J 



(Sn) ■ u t 



1 - 77 ) n ■ 1 



dS 



(4.17) 



bu t 



dx, 



+ P j T,dx = j 
n n 

where (3 — const > 0. Hence, since a > 0, it follows that if the external sources vanish 

6 = 0, . 9 = 0, 
and if the boundary conditions on S imply that 



(4.18) 



(Sn) ■ u t = 0, 



1- 77 )n q = 0, 
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then 

d f ( 1 



Q 

This provides the Lyapunov functional fn( e +^\ u t\ 2 — 0v) dx, which is nonincreas- 
ing on solutions paths. Let us notice that the boundary conditions (1.3) ensure (4-18). 
The identity (4-17) has been used in deriving energy estimates in Lemma 4.2. 

Our goal now is to derive further regularity properties from the energy estimates 
(4.10). To this purpose we use the regularity results for parabolic systems in Sobolev 
space with a mixed norm, stated in Lemmas 3.4 and 3.5. 

Let us consider the viscoelasticity system (1.1) with boundary and initial condi- 
tions (1.3)i, (l-4)i, expresses in the form: 

u tt -Qu t = V- [A 2 (e-9a)] + b in Q T , 
(4.19) u = on S T , 

u\ t=0 =u , u t \ t=0 = u 1 in 0, 

where Q = Q 1 is the viscosity operator (1.7). 

Applying Lemma 3.4 (i), (ii) to system (4.19) we deduce, in view of the bound- 
edness of A 2 and a, the following 

Corollary 4.5. Let us assume that 

Ul eB 2 p - 2 ^(n), fe€i M (ft T ), p, (re (l,oo), 

and if 2 — 2/a — l/p > then the compatibility condition uo\s = holds. 

(i) If £ g Lp^(n T ) and 9 g L Pta (fl T ), p, a g (l,oo), then the solution u to 
problem (4-19) satisfies 

Ikt'lU^Cfi') < c IK'H w i.V2 fnn < c(||e|| £p CT (nt) + ||0||i Pia (fit) 

(4.20) 

+ l|6||ip, < ,(n*) + ||«l|| B 2-V<r (n) ) 

for t g (0,T] ; with a constant c depending on il, S, T, p and a. 

(ii) If We g -L PiCT (f2 T ) and V0 g L p . a (Q T ), p,cr g (l,oo), then the solution u to 
(4-19) satisfies 

(4 21) " £ *'" M ^ /2 ( nt ) " c||w *'H^( fit ) " c (H Ve lk.*(n*) + l|V^||i Pi<7 (n*) 

+l|6|k„,(n*) + ||«i|| B ;-Va ( n)) 

/or t g (0,T] ; wit/i a constant c depending on Q, S, T, p and a. 
Using (4.10) in (4.20) for p = 2 and a arbitrary finite we have 

ll £ t'lli2, CT (nt) <c(A) + ||«i|| B 2-2/., f2) + ||6|| £2CT(n t)) 
(4.22) 2 ' ff v ; 

= cAi(cr), o-g(l,oo), t<T, 

where A is defined in (4.10). 

For further purposes (see the proof of Lemma 4.7) we prepare now some inequal- 
ities between the norms of u and 9. 



18 



I. Pawkow and W. M. Zajczkowski 



Let us consider viscoelasticity system (4.19) rewritten in the form 

u tt - Qu t - Qu = V • [{A 2 - Ai)e - 9A 2 a] + b in tt T , 

(4.23) u = on S T , 
u\ t =o = u , « t |t=o=«i in 0, 

where Q = Q 1 . 
Lemma 4.6. 

(i) Let u € iTj(n), Mi € L 2 (f2), b e £ 2 (^ T ) and e £ 2 (^ T ). Tften a 
solution u to system (4-23) satisfies 

Hwt'lk.ooCn*) + IIQ 1/2m II£2,oo(o') + IIQ^'lU^n*) 

(4.24) < c (t)(||(9|| L2(nt) + ||Q 1/2 m ||x 2 (o) + ||«i|U 2 (n) 

+ ||6|U 2(nt) ) for t < T, 

with a constant c(t) exponentially depending on t. 

(ii) Let u € H 2 (n) n Hl(Q), Ul e Ho( fi )» b e V6> e L 2 (fi). 
TTien a solution u to system (4-23) satisfies 

IIQ 1/2M t'll£ 2 ,oo(o*) + IIQ«IU 2l00 (n*) + IIQ«i'IU 2 (fi') 

(4.25) < c(t)(||V^|U 2 (o t ) + ||Q 1/2 wi||i 2( o) + ||Q«o|U 2( n) 

+ ll&lk(n*)) ^ t<T, 

wii/i a constant c(t) exponentially depending on t. 

Proof, (i) Multiplying (4.23) i by u t , integrating over £1 and using the boundary 
condition (4.23) 2 gives 



(4.26) 



lj t j(\ut\ 2 + \Q 1/2 u\ 2 )dx + J \QV 2 u t fdx 
J [(Ax - A 2 )e + 9A 2 a] ■ e t dx + J b ■ u t dx. 



n 

Using the estimate 

j {(A! - A 2 )e + 6A 2 a] ■ e t dx <5 X J \Q 1/2 u t \ 2 dx + c(l/<fi) J [\e\ 2 + 6 2 )dx, 
q n n 

which results on account of the Young inequality and (1.16), we conclude that 

j t j{\u t \ 2 + \Q 1 l 2 u\ 2 )dx + j \Q^ 2 u t \ 2 dx 

(4.27) 

< c 



;J(6 2 + \b\ 2 )dx + Cl J (\u t \ 2 + \Q 1/2 u\ 2 )dx, 



where we distinguished the constant c\. 

Hence, omitting the last integral on the left-hand side, it follows that 



d_ 

It 



J {\u t \ 2 + \Q 1/2 u\ 2 )dxe~ Clt < ce' Clt J {9 2 + \b\ 2 )dx, 
n n 
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(4.28) 



J {\u t \ 2 + \Q 1/2 u\ 2 )dx < ce Clt J (9 2 + \b\ 2 )dxdt' 
+ e c ^ J(\ Ul \ 2 + \Q^ 2 u \ 2 )dx. 



Now, using (4.28) in (4.27) and again integrating the result with respect to t' G (0, t) 
leads to 

J (\u t \ 2 + \Q 1 ' 2 u\ 2 )dx + J \Q x l 2 u t \ 2 dxdt' 

< c{te Clt + 1) J (d 2 + \b\ 2 )dxdt' + {te Clt + 1) J {\ Ul \ 2 + \Q 1/2 u \ 2 )dx, 
which proves (4.24). 

(ii) Multiplying (4.23)i by Qu t , integrating over fl and integrating by parts yields 



ld_ 

2 dt 



(4.29) 



J(\Q 1/2 u t \ 2 + \Qu\ 2 )dx + j \Qu t \ 2 dx 
= J (V • [(Ai - A 2 )e + 9A 2 a}) ■ Qu t dx - j b ■ Qu t dx = TZ. 



In view of the boundedncss of A\, A 2 and a, 

K<6 2 J \Qu t \ 2 dx + c{l/6 2 ) J (|Ve| 2 + \V6\ 2 + \b\ 2 )dx. 



q n 

Hence, choosing S 2 suffciently small and recalling the ellipticity of the operator Q, we 
get 

d_ 

dt 

(4.30) 



■ J(\Q 1/2 u t \ 2 + \Qu\ 2 )dx + J \Qu t \ 2 dx 



< c j (\\79\ 2 + \b\ 2 )dx + c 2 j (\Q 1/2 u t \ 2 + \Qu\ 2 )dx, 

n n 

where we distinguished the constant c 2 . Omitting the last integral on the left-hand 
side the latter inequality leads to 



d_ 

lit 



(\Q 1/2 u t \ 2 + \Qu\ 2 )dxe 



-c 2 t 



<ce~ C2t J (\V0\ 2 + \b\ 2 )dx, 



which after integrating with respect to t' £ (0,t), leads to 

J (\Q 1/2 u t \ 2 + \Qu\ 2 )dx < ce C2t J(\W9\ 2 + \b\ 2 )dxdt' 

J (\Q 1/2 Ul \ 2 + \Qu n \ 2 )dx. 



(4.31) 



+ e 
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Inserting (4.31) into (4.30) and again integrating the result with respect to t' G (0,i) 
gives 



J(\Q 1/2 u t \ 2 + \Qu\ 2 )dx + j \Qu v \ z dxdt' 

n 

< c(te C2t + 1) J (|W| 2 + \b\ 2 )dxdt' + (te C2t + 1) J (|Q 1/2 Mi| 2 + |Qu | 2 )^ 

q* n 

which proves (4.25). □ 

^From (4.24) we conclude that 

IK'IU 2 ,oo(n*) + l|w|| Loo (o,t;Jfi(n)) + \\Ut>\\L 2 (0,t;Hi(n)) 

< c(*)(||0||l 2( q*) + ||«o|Ui ( n) + ll«i|U 2 (n) + ll&IU 2 (n«)), * < T. 
Similarly, from (4.25) it follows that 

(4 33) ll^t'lliooCO.t;^ 1 ^)) + W u t'\\L 2 (o,t-H 2 (n)) + INL^O,*;*? 2 ^)) 

< c(t)(\\W9\\ L2{nt) + \\u Q \\ H 2 in) + \\ Ul \\ HHn) + ||6|U 2( Q t) ), t < T. 

Hence, by the definition of e, 

^ ^ \\St'\\v 2 (iV) = lkt'IU 2 ,oo(n*) + IkHU^o,*;*? 1 ^)) 

< c(t)(\\V9\\ L2{nt) + \\u a \\ H 2 {Q) + ||txi|| ff i ( o) + ||6||i 2 (a*)), 



where t < T. With the help of this inequality we prove 

,10 



Lemma 4.7. Assume that u G H 2 (Q), Ui G B 2 2 ^ 5 (fl), 9 G L 3 (Q), b G 



£ 2 ,io(^), ffG L 2 ,i(n*), * < T - 
T/ien 

(4.35) ll^i)IUs(O) + l|e|| ia (o,tiffi(n)) + IMIv 2 (o<) < cA 2 , t < T, 

w/iere 

A 2 = ip(Ao, Ai(10), ||uo||jj2 (n ), ||«i||jri(n), ||^o|U 3 (n), ||b|U 2 (n*), IMU 2jl (n*)) 
' < ¥>(H«o||ff2(fi), ||Mi|| B 2-i/5 (fJ) , ll^olkato), l|b|U 2 , 10 (n*)i IMU 2 ,i(n*)) = A 3, 

wii/i j4i(-) defined in (4-22) and Aq in (4-10). 

Proof. Multiplying (1.2) by 9 and integrating over f2* we get 



l|0(t)lli, ( n) + / \W)\\\ 2{n) dt> 

(4.37) 



< c 



9 2 \e v \dxdt' + c J J 9\e t >\ 2 dxdt' + c J J 9\g\dxdt' + ||0 O |||. 



(il)- 



on on on 

With the use of the Holder inequality the first term on the right-hand side of (4.37) 
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is estimated by 

t 

J ||0|U,(n)ll*lk(n)ll e *'IU a (n)dt' 



<sup||0(t)|U 3(n) J \W)\\ L6{ n)\\e v {t')\\ L2{n) dt' 
o 

< sup ||«(t)IUs(n)II^IUe. a (n*)He*'IU 3 (n*) 

<*! sup ||fl(t)||i, (n) + c(i/<y 1 )||e||^ (nt) || et ,||^ nt) = h. 

Using (4.22) for a = 2 yields 

h < 5isup||(9(t)||i s(n) + c(l/6 1 )\\e\\%\ 0<t;H1{n)) A 3 1 /2 (2). 

The second term on the right-hand side of (4.37) is estimated by 
II^IU e (n)||et'IUs(n)l|et'IU 2 (n)<ft' 



< ll^llL 6 ,2(f2 t )ll £ t'lli3,a 1 (0*)l|£t'll£2, CT1 (f2 t ) 
< 



"1 \ 

S2\\0\\l 2 ^ t , HHn)) + c(l/J2)||et'|li s ,. 1 (n*)^i^i). 



where (4.22)was used with a = ai, — + — = 2, Si > 2 but close to 2, because o\ can 
be an arbitrary positive finite number. 
Now we examine the integral 

si/3 \ l/si / \ \ l/si 



00 

(| lkf(*')l£ 1 (n) /8 *') 1/ ' lAl (/ IMOII£$ 3 *' 



" I / ik/'(/')iir <n> \ I y IMOIIS^'J = ^ 

where 1/Ai + 1/A 2 = 1. Setting siAi = 3, S1A2 = ^157, we obtain 

y imoiiW*'] (j n^(t')iii^)^ 
v 

Now, using (4.22) with cr = gives 

^<c(| ||^(i')ll^(o)^) 1/3 ^i /3 (^ 


for any si € (2, 3). 
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Finally, the third term on the right-hand side of (4.37) is estimated by 

c l|0|U2,oo(fi*)ll5lU 2il (n*) < cA \\g\\ L2 l{nt) . 

Inserting the above estimates into (4.37) and assuming that Si, S 2 arc sufficiently 
small, we arrive at 

IW*)lli,(n) + \\nUo,f,m (a) ) < P\\l 2m + cAl / \2)\\6\\l% ;Him 

(4 ' 38) +^ 1 (^)4 /3 (^)lkH| 4 / 2 (o,^(n ) ) +cA \\g\\ L2 , im + ||*o||i 8 (n). 

where — + — = ±, 2 < s 1 < 3. 

S± (7 1 Z 1 

Let us choose s\ = |. Then = 5 and o\ = 10. Since A\(a) is an increasing 

function of a it follows from (4.38) that 

L 3 (0) + \\Q\\L 2 (0,t;H 1 (n)) - M^lliaCn*) 



, + ^HL(0,t;Hi(n)) < ll^ll 2 

+ C (^(10)+< 3 (10)||e t ,|| 4 L / i t;ffl(f2)) + A)||.9||l 2>1 (q*)) + l|0o||! 3 (n), 

where we used the Young inequality in the second term on the right-hand side of 
(4.38). 

By virtue of (4.34) and (4.10) we obtain from (4.39) the inequality 



lkt'lk 2 (Q<) < c(t)(A + Al(10) + Ai(10) + Al /2 \\g\\^ Ant) 
+ \\0o\\ 



Employing (4.40) in (4.39) yields (4.35). This completes the proof. □ 
Using (4.35) in (4.33) implies 

, . lkt'||.L oo (0,t;.H-i(n)) + \\ u t'\\L 2 {0.t:H 2 (n)) + \\ u \\L^(a,t;H 2 {Q.)) 

< c(t)(A 2 + ||uo|U2 ( n) + 1 1 "ilk 1 (52) + l|b||x 2 (n*)) < c(t)A 3 , 

with A 3 defined in (4.36). 

^From what it has already been proved we deduce 

Lemma 4.8. Assume that u a € H 2 (fl), m e B 2 2 ~^(Q) n B 2 2 ~ 2 ' a {Vt) , 8 e 

H 1 (n),beL 2 . 10 {n t )nL 2ia (n t ), g eL 2 (n t ),a>4. 

Then 

(4.42) INk(n<) + l|V0|| Loo( o, t ;L 2 (o)) < ¥>(Ai(ct)), a > 4, 

w/iere 



Ai(ct) = ||«o|U2 (n) + ||ui|| B2 -i/5 (n) + ||«i|| B 2- a/ a (n) 

+ l|fo||/fi(n) + 11^11x2,10(0') + l|b||L a ,a(n*) + IMU 2 (fi<)- 
Proof. Multiplying (1.2) by 9 t and integrating over f2 gives 

c v J 0%dx+~ J \V0\ 2 dx 
n n 



(4.44) 

< c 



n 



0\e t \\O t \dx + c J \e t \ 2 \9 t \dx + J \g\\0 t \dx. 



3 (Q)- 
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Applying the Holder and the Young inequalities the first term on the right-hand side 
of (4.44) is estimated by 

eefd^j ' (^J 9\e t \ 2 dx^ ' J 09 2 dx + c{l/S 1 ) J 9\e t \ 2 dx, 

n n o n 

where on account of (4.35) the second integral is bounded by 

c (y e^dx^j 1 | £t i 3 ^ 1 <cA 2 \\£ t \\i 3 

n q 

The second term on the right-hand side of (4.44) can be estimated by 

c||0t|U 2 (n)IMk(fi)||et|U 3 (fi) < *2||^||i 2(n) +c(l/<y 2 )||e*||i 
Finally, the third term on the right-hand side of (4.44) is bounded by 

5 3 \\0t\\l 2m +c(l/5 3 )\\g\\l 2{n) . 

Employing the above estimates in (4.44), assuming Si, i — 1,2,3, sufficiently small, 
recalling that 9 > 0» > 0, and integrating the result with respect to time, we conclude 
on account of (4.41) that 

ll^ll! 2( n t ) + l|V^)||| 2(n) <cA 2 ^ 

' ' +C^||£ t '|lL(0, t ;£ 3 (0))+c||.9llL( f2 t)+c||Vc?o||l 2(0) . 

In view (4.41) and (4.45), applying Corollary 4.5 to problem (4.19) we conclude 

that 

IK Wwl-^v) < < A 3 + A /2a 3 + Mid' lUooCctjisCn)) 
(4-46) + ||ff|| L2(n *) + ll&IU 2 , CT (n*) + \\0 o \\m(.n) + ll«i||^-v- (n) ) 

< cA 3 ||e t /|| Loo(0it;i3(n)) +cip(A 4 (a)), 

where a is an arbitrary finite number. 

By the definition of the tensor e, inequality (4.46) implies 

(4-47) ||e t ' \\ w i,i/2 {nt) < cA 3 \\et>\\ Loo{0tt . L3m + cA 4 (cr). 

In view of the interpolation inequality 

(4.48) ||et'|U oo (0,t;£s(n)) < S\\e v \\ w ^i/2 (ttt) + c(l/<5)lkHk,„(fi<)> 

which holds for a > 4, it follows from (4.47) and (4.35) that 
(4-49) \\e t > 11^1.1/2(0.) < <fi(A 3 , A 4 (a)), u > 4. 

Hence, 

(4-50) \\ et , \\ Loo (o,t;L 3m < <p(A3, A A (a)), a > 4. 
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Applying (4.50) in (4.45) gives (4.42). This completes the proof. □ 

Let us note that since t is finite, estimate (4.49) in conjunction with the Holder 
inequality implies that 



(4.51) 



\ £t 'Ww 1 2 - 1 / 2 (nt) - p(*>^3,Ai(ct)), 



where a > 4 and cro > 1. 

Similarly, by (4.46) and (4.50), it follows that 



(4.52) 



u t'Ww 2 2: l o (nt) ^ v(*> A 3,Ai{a)), 



where a > 4 and ctq > 1. 

Lemma 4.9. Assume that 9 e Lp i0 -(fi*) } b e Zp i0 -(fi*), u € Wp(fi), iti € 

B2- 2 /'(n), pe (i,oo),ae(l,oo). 
T/ien 



(4.53) 

where 



\ u t'Wwiy 2 (nt) - c (*)[l|0|U P ,„(n*) + 4s(p,<r)], 



A 5 (p,cr) = ||« ||wj(n) + ll u illBj- 2 /-(n) + H b IU P ,-(n*)' 



and i/ie constant c(t) depends exponentially on t. 

Proof. Let us consider system (4.19) and apply the inequality (4.20). Representing 
e by 



(4.54) 



r(t) = J £ V {t')dt' + e(0), 



and using the generalized Minkowski inequality, we obtain 



t t' 



e t »(t")dt' 



o o 



£ P (H) 



dt' + c| ||e(« )||2 p(o) dt' 



t t' 



< c 





t t' 



e t „(t")\\ Lp{Q) dt") dt' + ct\\e(u )\\% 



(fi) 



< c 

' 

Consequently, denoting 



/ (/ ll^^Oll^cn)^) (* / )"- 1 d* / + ct|| e («o)||£ p{ n). 



a(t) = ||et'|l£ p ,„(n*) = J \\ £ t>(t' j\\ Lp{n )"» . 
o 

we deduce from (4.20) the inequality 

t 



a(t) < J a(t')a(t')dt' + A(t), 
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where; 



a(t) = ct a -\ 

Mt) = c(\\6\\l pA nt) + \\b\\l pM nt) + *l|e(«o)ll^ ( n) + ll«ill^-v. 
Hence, by the Gronwall inequality, it follows that 



Z 



< A(t)(l + ai(t)e ai ^), ai(t) = ta(t). 



Thus, 



(4.55) 



+ H u i|lB-^ (n) + l|b|li t) ,.(n*)) 



= c 



with c(f) = c(l + a 1 (t)e ai ^). 

Using (4.55) in (4.54) yields the analogous bound on a (Q*)- Consequently, on 

account of (4.20) the corresponding bound on ||Mt'|lw 1,1/2 (fit) f°ll° ws as well. The 
proof is completed. □ 

On the basis of Lemma 4.9 we prove now 

Lemma 4.10. Assume that u G Wl(ty, Ui G B 2 r ~ 2/r (9), 9 G L r (Q), b G 

L r (ft T ), 3 G M nT )» r G (l,oo). 
Then 



(4.56) 
where 
(4.57) 



li,,,^*) < 4s (r, »",*), r<oo, 



A 6 (r,r,t) = c(t)(||0 o ||L r (n) + </fA a (r,r) + \\g\\ LASi *) + 1), 



and As(p, a) is defined in (4-53). 

Proof. Multiplying (1.2) by 9 r , r > 1, and integrating over f2 gives 



ir+2 



4fcr 



dx 



(4.58) 



c v d 

H2dt rlli ^) T (r + 1) 2 

= ~y 6» r+1 (A 2 a) -e t dx + ^ 6> r (A l£t ) • e t dx + ^ r 5 dx. 
n 12 n 



Hence, after integrating with respect to time, 



r+2 



Akr 



r + 2 n»wii Vl (n) ' (r + 1)2 



r+2 

V6»— 



dxdt' 



(4.59) 



S2' 



6> r+1 |e t ,|dxdt' + c y e r |e t /| 2 da:<ft' + J 6 r gdxdt' + Jjy o ||£J a(n) . 
n* n« n* 



26 



I. Pawkow and W. M. Zajczkowski 



Here let us recall Lemma 4.9 which provides the inequality 
(4-60) \\et'\\L A n*)<c(t)(\\e\\ Lr{nt) +A 5 (r,r)) 
for r G (1, oo). 

On account of (4.60) the second integral on the right-hand side of (4.59) is esti- 
mated as follows 



J ^|£t'| 2 ^'<||^||2 r+2(at) || e Hll r+2 (a*) 



(4 ' 61) < c(t)\\9\\l r+2{nt) (\\ef Lr+2m + A 2 5 (r + 2, r + 2)) 

<^)(ll^llt+ 2 (o t )+^ +2 (^ + 2^ + 2)). 

Now, using (4.61) we estimate the first integral on the right-hand side of (4.59) by 



/ 



e r+1 \e t ,\dxdt' = J e r 2 +1 9^\e t '\dxdt' 



< \\^ +1 \\ L2m W^'WLw - ii e iiZ" +2( ^)ii^i^i 2 ii^) 

< <m\T +2 (n^T +2 m + A ?^( r + 2 ^+ 2 )) 

<c(t)(l|e|r L t 2 2( ^ ) +^ +2 (r + 2,r + 2)). 

Finally, with the use of Holder's and Young's inequalities the third integral on the 
right-hand side of (4.59) is estimated by 



e r gdxdt> < \\e\\ r Lr+2{nt) \\g\\ L(r+2)M nn < 40\\Zl ( nn +«\\g\\?* am + !)■ 

Inserting the above estimates into (4.59) leads to 

t 

ll*(*)C 2( n) < ^ C W(/ ll*«^ ( n)^ + ^V + 2,r + 2) 

o 

+Ni2: 2 2( o t )+i)+n^ic ( o). 

Hence, by the Gronwall inequality, we conclude that 

WmilZlm < [\\Oo\\Z 2 +2m + (r + 2)c(t)(Al + \r + 2,r + 2) 
+ \\9\\Zl 2( n t) + l)}eMc(t)(r + 2)) 
<A r + 2 (r + 2,r + 2,t) 

for t e (0,T), with A 6 defined by (4.57). This gives (4.56). The proof is completed. □ 
Corollary 4.11. Taking into account that \fr is bounded let us define 

Mr,r,t) = c(t)(\\u \\ w i {n) + \\ui\\ B 2-2/r {n) + \\e \\L r (n) 
+ ||6|U r , r (n*) + \\g\\L r (nt))- 
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Then, according to (4-56) and the definition of A 5 (r,r) in Lemma 4-9, 
(4-63) l|0|k,oo(n*) <A 7 (r,r,t), re(l,oo), t <T. 

Moreover, by (4-53) and (4-63), 

lkHk,„(fi') < c\\iH>\\ w w m 

< c(t)A 7 (r,r,t), {r,a) € (l,oo). 

Let us consider now the elliptic problem resulting from (1.2) and (1.3)2: 
- kA6 = -c v eB t - 6{A 2 a) ■ e t + (Aie t ) • e t + 9 in fl, 

(4.65) 

v ; n-V6» = on S. 

We have 

Lemma 4.12. Asume that u G Jf 2 (fi) n (Q) n JS 2 ^ 27 ^), 

mi g B 9 2 %(n)nB^^(Q), e g fl^fi) ni^(fl), 

b g i 2 ao(^ T ) n V(n T )n%(fi r ), s e £ 2 (^ T ) n £_2s_(ft T ). 

Then 

(4.66) IMIw^di*) ^ ^(s.M)), * < T, 

w/iere 

A s (s,a,t) = c(t)(\\u \\ H 2 (n) + \\u \\ w i 2s (n) 



(4.67) 



"-«2,10( Si ) 2 , a ( Si J "b 2 J 2s (Q) 



2-s ' 2- 



+ ll^ollffi(n) + l|0olU^_ ( n) + ||&lk, 10 (ii«) + ll&lk,.(n<) 
+ \\ b h 2s (n*) + IMk(n*) + II.9IU 2s (a*)), 



and 1 < s < 2, ct > 4. 

Proo/. On account of (4.42), (4.63) and (4.64) it follows from (4.65)i that 

\\A6\\ L3{nt) < c(\\6\\ L ^_ {nt) \\6 t \\ L2{Qt) 

+ l|0|k s (n*)lk*'lk(n*) + ||et'HL(n*) + IMk(n*)) 
(468) < c(t) U( * ^ )t Wa)) 



2 - s 2 - s 

,(Q«) 



+ A?(2 s ,2 s ,t) + ||« 7 || i 



where s < 2 is close to 2, cr > 4 and t <T. 
In view of the estimate 

(4-69) ||0||w?(i>) < c(||A0|k (n) + ||0|k ( n)), 

which holds for the Neumann problem (4.65), we deduce from (4.68) that 



(4.70) 



■|k(o,t ; w?(n)) < cPU^n*) + c{t) (a 7 QH£_ tj <p(M<r)) 

+ A 2 7 (2s,2s,t) + \\g\\ L3{nt) y 
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where er > 4, s < 2 close to 2, t < T. 

Now, combining (4.70) with (4.42) and using (4.10) we conclude that 
\\0\\ w ^ m < c(t)(A Q + viA^a))) 

(4.71) +C WK^^'^ (i4(ff)) 



+ i4?(2s,2*,t) + ||0|| i<( ni)^ < y{A s (s,a,t)), 



where the latter inequality follows from the definitions of the quantities A (see (4.10)), 
At(cr) (see (4.43)), A 7 (r,r,t) (sec (4.62)) and Ag(s,a,t) (see (4.67)). This proves the 
lemma. □ 

Corollary 4.13. In view of the imbedding 
(4.72) VW^(n T ) C L p ,, s (n T ), ae(l,2), 

where s < p' < 3 ^-i > an< ^ ^ e Holder inequality with respect to the integral over VL in 
the case 1 < p' < s, we have 

(4-73) ||Vw|| V s( n r ) < c||u|| w? ,i (nT) 

for any v G W^(n T ), K p' < jj^. 

Applying the imbedding inequality (4-73) to (4-66) we deduce that 
(4-74) ||W|| Vs(nt) <c(t)<p(A 8 (s,o-,t)) 

where 1 < p' < , s G (1,2), a > 4, and As(s,a, t) is defined in (4-67). 

We shall use (4.74) to get more regularity estimates on e. To this purpose we 
return to the viscoelasticity system (4.19) and prove the following result analogous to 
Lemma 4.9. 

Lemma 4.14. Assume that V6» G L p . a (tt T ), b G £ p>CT (fF), 
Then 

W^'Wwi-v 2 ^) ^ c IML^(fi<) 

(4.75) < c(t)(||Vfl|U Pitr(n *) + ll&IU P ,„(n*) + ll«i|| B ;-v* (n) ) 

= c(t)(\\ve\\ LpAnt) + A 9 ( P ,cT)), 

where t < T , p,a G (1, oo), and 

A 9 (p,(j) = ||«i|| B 2-2/. (n) + H&ll^cn*)- 

Proof. Let us consider system (4.19). By inequality (4.21) we have 
||Ve t '||x Pi<T (nt) < \\u t >\\ 

(4.76) 

+ HHI^-v-^) + ll b IU P ,.(n')) 

for t <T. We use now the formula 

t 

(4.77) Ve(t) = y Ve v {t')dt' + Ve(0), Ve(0) = Ve(« ), 

o 
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and repeat the proof of Lemma 4.9 with Ve, V0 in place of e, 0. In result we conclude 
that 

(4.78) ||VeHk,„(Q<) < c(t)(\\Ve\\ Lpa{nt) + ||ui|| B; - V . (n) + ll&lk,^))- 

Using (4.78) in (4.77) implies the analogous bound on || Ve|| ip a (Qt) and then, by 
(4.76), on Hwt'llw^'iCfi') as wel1 ' D 

Corollary 4.15. Applying estimate (4-74) in (4-75) yields 

(4.79) ||e t '|| w W2 (nt) <c(t)(<p(A 8 (s,a,t))+A 9 (p',s)) 

for Kp'< 37^, s G (1, 2), a > 4, t < T. 

Corollary 4.16. Let us consider the imbedding 

(4.80) W 1 // 2 ^) d L^ 2 {n T ), s €(1,2), 
which holds true provided p' > ~r . This condition together with 

s 

p' < -j^Y implies that 

5 

-< S <2. 
Consequently, it follows from (4-79) that 

(4.81) \\e v |U 0O , a( n*) < c(t)(<p(A 8 (s, <j, t)) + A g {p' , *)) = A 10 (s, a,p', t) 

for ^ <p> <^L-, se(|,2), (7>4. 

Estimate (4.81) plays the key role in getting Loo(ri T )-norm bound for 9. 
Lemma 4.17. Let the assumptions of Lemma 4-12 be satisfied. Moreover, let 
tti G B 2 p 7 2 J s {n), d G Looifi), b G V,4^ T ), .9 e £oo,i(^ T ), w/iere 

A „\ 3 ,3 
•^(3,2), _ 1 <p'< 1 _ 1 , <7>4. 

TTien 

(4g2) ll^llLcCn*) <c(0[Aio(s,a,p',t)+Af (s, £ 7,p',t) 

+ ll»olk 00 (n) + IMIwfn*)] = c(iL4 n (s,<7,f/,i), 

where Aio(s,a,p' ,t) is defined in (4-81). 

Proof. Let us consider once more the identity (4.58). Regarding (4.81), we have 

(4.83) 7T2dt" 6, "^+ 2 (°) - c ll 6, llit!i(^)ll £t ll £ ~(") 

+ c||^|IL(n)ll £ *HL(n) + PWl^hWL^ny 
Taking into account that > (9* > 0, and using the inequality 
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and 



< 110117 < 



|r+l 



lL P (n) ^ l»l' T -|l"lli P+a (n) ^ IIMwn) r+a(n) 



where constant c is independent of r, we deduce from (4.83) that 

MZUn^ML^m < c\\e\\Zl 2{n) (\\s t \\ Loo{ ^ + \\s t \\ 2 Looin) + \\g\\ Leo{a) ). 

Hence, 

^|| 0\\ Lr+2 (Q) < c(||£t|| ioo (n) + ||£t|lL(o) + HfflUooCn)), 

which in conjunction with (4.81) leads to 

||0(*)IU r+2 (n) < ||^o|U r+2 (n) +c(||e t '|Uoc,i(n t ) + Ikt'llL^n*) 
(4-84) +\\g\\ Lool{nt) )<c(t)(A 10 (s,a,p',t) + A 2 w { S ,a, P ',t) 

+ l|0o||z, r+2 (fi) + II^IUoo, 

Now, letting r — >• oo in (4.84) we get the assertion. □ 

Corollary 4.18. Repeating the proof of Lemma 4-12 with the use of the upper 
bound (4-82) allows us to deduce that 

(4.85) \\B\\ w F m <<p{t,An{s,<T,j/,t)), 

where t<T,sE (f,2), ^ < p' < jj^j, a > 4. 
Corollary 4.19. Applying (4.82) in (4-53) yields 

(4-86) IkHU^n*) <c(t)[i4n(«,<r,j/,t)+4s(p,tr / )], 

where t <T, s e (§,2), 5^77 < p' < 3/s-i ' a > 4, p, a' G (1, 00), A 5 (p, <r) is defined 
by (4.53) and A n (s,a,p',t) by (4.82). 

The next step consists in obtaining a "better" estimate for 9 by means of the 
parabolic regularity result stated in Lemma 3.5. To apply this result to the quasilinear 
equation (1.2) we need to prove first that 9 is a continuous function on f2 T . We are 
able to prove more, namely the Holder continuity by means of the parabolic De Giorgi 
method in the same way as in [18, Lemma 6.1]. 

Since the above reference concerns much more general situation, we present here 
for reader's convenience a direct, simpler proof. 

Following [14, Chap. II. 7] we record the definition of the space 
B 2 (ti r ,M,'y,r,6,K), where fl T = Q x (0,T), ft C K", n G N, and M,-y,r,S, 
k are positive numbers. 
The function u G # 2 (fF , M, 7, r, 5, k) if: 

(i) «e^°(fl T ) :-C(0,T;L 2 (r!))ni 2 (0,T;H/ 2 1 (r!)), 

(ii) esssup^T |u| < M, 

(iii) the function w(x, t) = ±u(x,t) satisfies the inequalities 



max 

*0<*<*< 



>c +r \\(w - k) + \\ L2{Be _ aie{xo)) <\\{w- k) + (;t )\\ L2{Be{xo)) 
+ Tl^e)- 2 !! (w fc) + ||| 2(Q(eiT)) + ^d+^fc, o, r)] 
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and 

\\(w - k)+\\y 2{Qie _ aieiT _ a2T)) 

< l{[(°iQ)- 2 + (V2T)- I ]\\{w fc) + ||i a(0(BiT)) + ^ 1+K \k, e , r)}. 

Here the following notation is used: 

(w — k) + = maxfifi — k, 0} — the truncation of w, 

B e (xo) = {i £ SI : \x — xo\ < g} — a ball in fi, 

Q(£i,t) = S e (a;o) x (t ,t + r) = {(z,t) efi T : \x - x \ < g, 

t < t < t + t} — a cylinder in fl T , 

where £>, r are arbitrary positive numbers, a\, a 2 are arbitrary numbers from the 
interval (0, 1), and k is an arbitrary number such that 

ess sup w(x,t) — k < S. 
Q(g,r) 



Moreover, 



V 2 (n T ) = L^O, T; L 2 (fi)) n L 2 (0,T; W^fi)), 



t0+T 



fi(k,g,r)= J mc&s r / q Ak, e (t)dt, 

to 



where 



AkAt) = i x e ^efco) : w(x,t) > k}, 

and positive numbers q, r are linked by the relation 

Inn 
r + 2q ~ 4' 



with the admissible ranges 

2n 



qe 2 



€ [2, 00) for n > 3, 



n - 2_ 

ge(2,oo), re (2,oo) for n = 2, 
ge(2,oo], re [4, 00) for n = 1. 

Lemma 4.20. Lei ifte assumptions of Lemma 4-17 be satisfied, and there exist 
constants and M suc/i that 

9>9*>0, 
M=\\e\\ Loo{nT) <c(T)A n (s,a,p',T), 

IMImo-) <c(T)[i4 11 («,o-,j/,T)+A B (p,p)] ) 
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where p G (l,oo), An(s, a,p', T) is defined by (4-82) and A 5 (p,p) by (4-53). Besides, 
let k be a positive number such that 

(4.88) snp9 (x)<k 

o 

and 

M — k < 8 with some 6 > 0. 

Then 

(4.89) ^B 2 (H T ,M, 7 ,r,i,K) 
with 

10 / 2 

r = <z = Y' K H'3 

^c { t)(l + A 2 l{ s^,T )+ A 2 ^^) 
, 2 / 10 10 \ 

Proof. Note that the bound (4.87)i is ensured by Lemma 4.1, the bound (4.87)2 
by Lemma 4.17 and (4.87)3,4 by Corollaries 4.18 and 4.19. 

We check that 9 satisfies the conditions (i)-(iii) in the definition of Z?2(^ T , M, 7, r, 5, n). 

^From [1, Chap. 3, Sec. 10] we have the imbedding 
W 2 2 ' 1 (ft T ) C V r 2 1,0 (fi r ). Since 9 e VK 2 4 (0 T ) the condition (i) is satisfied. 

Condition (ii) is automatically satisfied on account of (4.87) 2 . 

Let us check that 9 satisfies the second inequality in condition (iii). By virtue of 
(4.87)i it sufficies to consider (iii) with w(x,t) = 9{x,t). 

Let Q(q,t) = B e (x ) x (to,*o + t) be an arbitrary cylinder in fi T , and Qix^) 
be a smooth function such that suppC(x,i) C Q(q,t) and C( x ^) = 1 f° r ( x ,t) G 
Q(q — aiQ, t — cr 2 r), where a\, cr 2 e (0, 1). Moreover, let 

^fc, e (*) - {a; e B e (a;o) : 0{x,t) > k}. 
Multiplying equation (1.2) by ( 2 (9 — k) + and integrating over Q, gives 

y y 9( 2 d t (9-k) 2 + dx + k J \V(9-k)+\ 2 ( 2 dx 

(4.90) Q ° 

+ 2fc / CO? - fc)+V(0 - fc)+ • V(dx = / GC 2 (0 - k) + dx, 
n n 

where, for simplicity, the right-hand side of (1.2) is denoted by 

G = -9{A 2 ol) ■ e t + (Aie t ) • s t + g, 

and to avoid the notational collision the letter k for heat conductivity is replaced by 
ko (for this proof only) . 



c 
~2 

(4.91) 
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Let us rearrange the first integral on the left-hand side of (4.90) to the form 
I J 9 C 2 d t (9 - k)ldx = %j t \w ~ k) 2 + C 2 dx 

- ^ J 6 t {6 - k)\C 2 dx -c v J 9(9 - k) 2 + (( t dx. 
n n 

Further, the term with 9 t in (4.91) is rearranged as follows 

- I / t (9 k) 2 + ( 2 dx = -5l[(0- k) 2 + d t (9 k) + C 2 dx 
n n 

(4.92) =-jJ d t (9-kf + ( 2 dx 

n 

= ~j t j{e - k)lc 2 dx + jJ(o- k)lcc t dx. 

n n 
Inserting (4.91) and (4.92) into (4.90) we obtain the identity 

J 9(9-k) 2 + ( 2 dx + k J \V(9- k) + \ 2 ( 2 dx 
n 

•- k )\c 2 dx- c ^ J(9-kf + a t dx 

•„ J 9(9 - k) 2 + (( t dx - 2k J ((9 - k) + V(6 - k)+ ■ \7(dx 
+ J G( 2 (9-k)+dx. 

Q 

Integrating (4.93) with respect to time, taking into account (4.87) i and the fact that 
by (4.88), (9q — k) + = 0, we conclude that 

^ J(9-k) 2 + ( 2 dx + k J \V(9-k)+\ 2 ( 2 dxdt' 
<cj (9-k)\C 2 dx + c j (9- k)\\C\ \C,t>\dxdt' 

J 9(9 - k) 2 + \(\\Ct>\dxdt' + c J (9 - k) + \V(9 - k) + \\(\\V(\dxdt' 

J \G\(9 - k)+( 2 dxdt' = Y, J i- 



c v d 
~2~dt 

a 

c v d 
" ~6di 

(4.93) 

< 

n 



(4.94) 



+ c 
n 



Since, by the assumption M — k < 6, it holds 

(4.95) (9 - k)l < 5(9 - k)\ with arbitrary 5 > 0, 
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the integral 1\ can be absorbed by the left-hand side of (4.94). Further, using (4.87) 2 , 
(4.95) and the bound |C| < 1, 



Next, by Young's inequality 



h + h < M J {9 - k)\\C t ,\dxdt' . 



h<^J |V(0 - k) + \ 2 ( 2 dxdt> + JLj(0- k) 2 + \V(\ 2 dxdt>, 

a* 

so the first integral on the right-hand side of the latter inequality is absorbed by the 
left-hand side of (4.94). In result, incorporating the above estimates into (4.94), we 
arrive at 

J(9- k) 2 + ( 2 dx + J \V(9-k)+\ 2 ( 2 dxdt' 
<c(M + l) j (9-k) 2 + {\VC\ 2 + \Ct>\)dxdt' + 



(4.96) 



where 



h = c 



J \G\(9-k)+( 2 dxdt'. 



By the definition of (, it holds 

M j(9-k) 2 + (\V(\ 2 + \( t ,\ 2 )dxdt' 

v 

M[(a ig )- 2 + (a 2 r)- 1 ] J (9 - k)\dxdt' '. 



(4.97) 

< 



Q(e,r) 



It remains to estimate the integral 1$. Recalling (4.87)2 again and applying 
Holder's inequality yields 

t +T 

2 dxdt' 



h = c J j \G\{6-k) + t 

to A k , e (t') 

l«+T -I t +T . 

(f f \ ' / / \ i / A 2 

/ / \G\ Xl dxdt'\ ( / mcasA k ^(t')dt' 

where 1/Ai + 1/A 2 = 1. To satisfy the conditions in the definition of the space 
^2(fi T , M, 7, r, 5, k) we set 

12 13 3 

— = -(1 + k) and r = q with — h — = -. 
A 2 r r 2r 4 
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Then 



A 2 = 



Ai = 



e l|oo 



for k € ( 0, - 



3(1 + k)' " l 2 -3k 

Consequently, 

I 5 <M\\G\\ Lxi{nT) ^ x *(k,Q,T). 
By virtue of (4.87)2 and (4.87)4, we obtain 

WGL^inn < c(\\0\\ LooinT) \\e t \\ L ^ {nT) + INI?.,, • 

+ HfllU j_(nT)) < c(T) [A 2 u (s, a,p', T) + A 2 ( 

2-3k \ \ 



2 - 3k 2 - 3k 



10 



10 



2- 3k 2- 3k 



+ IMU 5 (nn 

2-3k 



7i 



Hence, 
(4.98) 



h <M 7lM f( 1+K )(fc,g,T). 



In result, applying estimates (4.97) and (4.98) in (4.96) leads to 



\\(o-k) + \\ 2 V2iQ{ 



e-cie,T-cr 2 r)) 



ess sup ( [6 -k)\C 2 dx+ [ \\7(6-k)+\ 2 ( 2 dxdt 
te[o,T]J J 



< i{[(°iq)- 2 + {°2T)- l ]\\{e fc) + ||| 2(Q(eir)) + M « (1+K) (fc, e, r)}. 

This proves the second inequality in condition (iii) with the positive number 

(4.99) 7 = c (T)(Al 1 (s,a, p', T) + 1 + 7l ). 

The first inequality in (iii) follows by multiplying (1.2) by Q(0 — k) + , where C,q(x) 
is a smooth function such that supp(o(x) C B g (xo), Co(x) = 1 for x G B g - aig (xo), 
G\ G (0, 1), and next integrating over x (£ , to + r). Then, by repeating the presented 
above estimates we arrive at the following inequality in place of (4.96): 



J (9-k) 2 + C 2 dx + J \V(9 -k) + \ 2 Q dxdt 



(4.100) 



< c(M + 1) 



Q(e,r) 



(6(h) - ky+Qdx + j (6 - ky + \V( Q \ 2 dxdt 

B e (x ) Q(q,t) 

\G\(6 - k)+( 2 dxdt. 



Q{e,' 



The last two integrals on the right-hand side of (4.100) are estimated respectively by 
(4.97) with (ct2t) _1 = 0, and by (4.98). In result, (4.100) provides the first inequality 
in the condition (iii) with 7 defined in (4.99). The proof is complete. □ 
Corollary 4.21. By virtue of the imbedding (cf. [14, Thru II. 7.1]) 

B 2 (n T , M, 7, r, 5, k) C C a - a/2 (n T ), a G (0, 1), 
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2 - 3k 2 - 3k 



it follows from (4-89) that 

(4.101) 9 e C Q ' Q/2 (fi T ) 

wifft Holder's exponent a e (0, 1) depending on M, 7, r, <5 onrf k, where 

M = supfl < c(T)A u (s,a,p',T), 
7 (s, k, T) - c(T) (l + A 2 n (s, a,p', T) + A\(^ 

+ A5 ( v2 -^'2-3^J +ll5ll ^H' 
/5 „\ 3 ,3 A / 2 

S H3' 2 J'^A <P< 37^T' a>4 ' K H°'3 



In view of Holder's continuity of we can apply Lemma 3.5 and deduce the final 
estimates on 6 and u. 

PROPOSITION 4.22. Let the quantity 

Mi) = KUi^ 2/5 (n)rw;,(fi) + ll ni llBiV» (n)nB^"'(n) 

( 4 - 102 ) + ll^o||wi.(n) + IMIw^n*) + l|6|U ia (n«)n£ p ,, a ,(n«), 

p* > 3, e (1,00), 



be finite. 

Then the following a priori estimate 



(4.103) \\0\\ Loo{ n*) + \\e\\ c «.«/* {a >) + Ikf IU p ,, a<( n«) < <^(*)), * < T, 



is wo/zd. 

Proof. iFrom (4.82), (4.86) and (4.101) it follows that the estimates on ||0|| Loo ( n t), 
||et|U Pi „,(n*) and \\d\\ ca , a/ 2 (n t ) involve the parameters s e (f , 2) ,p' e (^§77' 37!=^ ' 
cr > 4, (p, a') G (l,oo) and k e (0, |). We select these parameters in such a way to 
express the above mentioned estimates in an explicit, sufficiently simple way. 
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Let us recall that the estimates depend on the quantities: 

A 8 (s,a,t) = c(t)[\\u \\ H 2 (n) + \\u \\ w i^ {n) + ||«i|| B 9/B o(n) 



+ «ib-v. 



«1 3s=2 



l|0o||ifi(n) 



(4.104) 



+ Pnh^n) + WbhzMW) + ll b IU 2 , CT (n<) + ||b||i^_(n*) 
+ IMU 2 (n*) + HfflU^cn*)], 

A 9 (p',s) = ||«l|| B 2-2/. ({1) + ll&Hi s(0 t), 

p' ,8 

A 10 (s, <j,p', t) = c(t)(ip(A 8 (s, a, t)) + A g (p', «)), 
A n (s,a,p',i) = c(t)(A w (s,a,p' ,t) + A 2 10 (s,a,p' ,t) 

+ II^oIIl^^) + hh^Aw))' 

j(s,a,p',K,t) = c^l + ^ 1 (s,o-,p',t) + A\(^ 



2 - 3k' 2 - 3k 



10 10 



2 - 3k 2- 3k 
5 



-4, 



2 - 3k' 2 - 3k 

+ H b H^^_("')' 
10 10 



2 - 3k' 2 - 3k 
+ ll&IU 10 (n), 



2-3k ' 2-3ft 



l"o||wi 10 (n) + ||«i|| sdga 

2-3« 13 10 10 

2-3« ' 2-3k 



A 5 ( P ,(t') = ||uo|| w i(n) + ||mi|| b 2-2/„' (0) + ||b||i P!(r ,(ot), 



where s G (|,2), p' G ^gTX, , a > 4, k G (0, §). 
To have the above quantities finite we need 

M e H 2 (n) n w\^(n) n w\o_(fi) n wi(fi), 



,2-2/ct 



2-2/ s. 



(4.105) 



2-3n'2-3n 2-3h'2-3k K ' 

b g L 2 ,io(fi*) n i 2 ,a(^*) n (n*) n Lp/, a (n*) n (o*) 
nL Pi<7 ,(n*), 

5 g L 2 (n*) n Z/_?s_ (n*) n Loo,i(^*) n £ 5 (O*), 



where s G (f ,2), p' G (^r, -^1 , cr > 4, k G (0, §). 

V S 8 J 



Let us set 
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Then £- s = 12, 2^ = 3, = 6, and (4.105) takes the form 

M g H 2 (ft) n w{ 2 (n) n w*(fi), 

m g B»/ 5 (n) n Bi^Vfi) n B& 6 (fi) n B*/ 6 2/7 (fi) n B^ /<T '(fi), 

(4.106) JoGioolfiJnff 1 ^), 

.9 e Loo,i(fi*) n Li 2 (n*) = £00,12(0*). 

Using the imbeddings 

^ 2/5 (S])cff 2 (!])n<(fl), 
^>)c<(fi)nBj /7 (fi), 

W^ft) C Loo (fiJniT^fi) for p>3, 

we replace (4.106) by 

«o e w? 2/B (fi) n w*(fi), mi g b^ 6 2 (Q) n sj-^fi), 

(4.107) o ew£(fi), p*>3, 

6 g Li 2 (n*) n £ p ,a'(fi*), .9 e ioo,i2(^ 4 )- 

The assumptions (4.107) ensure that the quantity A(t) in (4.102) is finite. Then 
estimates (4.82), (4.86) and (4.101) imply (4.103). Let us note that p is replaced by 
p'. This completes the proof. □ 

Lemma 4.23. Assume that the data satisfy (4-107), SeC 2 , and moreover, 

u e wl q (n), «i g B*-/j/f(fi) n B^(n), 

(4.108) fl g Bj- 2 /*^), 6 G L PlP0 (fi*) n L 29 , 290 (^), 
jel M0 (Sl'), 

wif/i any numbers p,po,q,qo G (l,oo) smc/i i/iai 

3 2 3 2 



< 1. 

9 9o P Po 



Then 



w, 



tUnt) < V (A(t),t)+cmi 2g , 2mm 



(4.109) + lluollw^n) + ll Ul ll B ^°(fi) + HfflU., w (n*) 

+ 1 |0 O 1 1 B 2-2/<, 0(n) ) = SiO?, 9o, 

and 

(4.110) lht|| w 2,i o(at) < c(t)Bi( g ,go,t)+c(||6|| Xppo(nt) + ||«i||^-v»o (n )). 

where A(t) is defined by (4-102), and t < T . 

We remark that, by Lemmas 3.4 and 3.5, the constant c in (4.109) and (4.110) 
depends on p, p , q, q and C 2 -norm of the boundary S. 
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Proof. In view of the Holder continuity of 9 along with its lower and upper 
boundedness we deduce from Lemma 3.5 and Proposition 4.22 that 

P\\ W 2 ; i {nt) < c(\\9(A 2 a) ■ ef\\ Lti9o(n *) + IK-Mf) ■ etflU ll40 ( nt ) 



+ ||ff||i liW (n«) + ll^o|| B 2 -2 



2/90 fQV 

(4.111) -"™ 1 ' „ 

< V(^(*),t)l|ef IU llW (n*) +c(||et'|lL 2 ,, 2 , (n«) 

+ ||5lU liW (n*) + II^oIIb^o^), t < T. 

^From (4.103) with the parameters (r,r ) = (p, a') selected in Proposition 4.22, we 
have 

(4.H2) l|et'IU r , ro( n*) <<p{A{t),t), 

where (r, r ) is equal (2g, 2<7o)- 
Using (4.112) in (4.111) we conclude (4.109). 
Let us consider now the imbedding 

v^,i(fi T )ci^(fi T ), 

which holds under the condition 

(4.H3) - + ---,-\<h q'>q, %>Qo- 

q qo q' <? 

Then it follows from (4.75) and (4.112) that 

\\ u t'\\wli (nt) < c(t)(\\Ve\\ LppQ{nt) + ||6|U PiP0 (n«) + ||ui|| B j-2/po (n) ) 
(4.114) < c(t)(||V6»|U 4li4 , (nt) + ||&|| W n«) + \\ui\\ Bp -v P0{n) ) 



< c w(\\ \\wll o (nt) + l|6|U P , P0 (n*) + ll M ill B ^o (a) ), 

where p < q', po < q' , and q',q' satisfy (4.113). Finally, we set p = q', p = q' Q . This 
establishes (4.110) and thereby completes the proof. □ 

Remark 4.24. Since u = f*u t >dt' + u , u t e W 2 '* o (Q T ), p,p £ (l,oo), and 
Mo G Bl- 2 J ao (Cl), it follows that u e W 2 p ^ o (fl T ) where a G [l,oo). 

COROLLARY 4.25. Comparing assumptions (4-107) with p = p' and (4-108) we 
deduce that 

Boe< /6 (fl)nwJ,(n)n^(n), 
(4 ' 115) u! g B\l%(n) n B 2 - 2 J°\n) n n 

Setting p' = 2q, a' = 2qo, q = qo = 6, p = po = 12, we get 

u g w 2 2/5 (n), ui g b^ 6 2 (0), 

and 

O g Wp 1 . (O) n B 6 v 6 3 (fi) = B 5 6 ^(n) for 3 < P , < 6, 
6 g £12^') n L PlP0 ( fit ) n L^^n*) = L 12 (n l ), 
g g Loo,i2(^*) n = Loo,i2(fi*). 
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Introducing the quantity 

B{t) = IKHw^ 2/5 (n)rw; 2 (n) + IMIsJv^n) 

+ l|b|U 12 (n*) + l|0o|| B 5/3 (n) + \\g\\ LooA2 (W), 
we conclude in place of (4-109) and (4-110) the estimate 
(4.H6) WMwtfw + PWw^tt*) < v{B{t)). 



5. Existence (proof of Theorem A). To prove the existence of solutions to 
problem (1.1)-(1.4) we use the following method of successive approximations: 

u n tt +1 - V • (A 1 e{u n t + 1 )) = V • [A 2 e(u n ) - {A 2 ct)6 n \ + b in Vl T , 
c v 9 n 9? +1 - kA6 n+1 = 

-6/ n (A 2 a)-e«) + A 1 £«)-£«)+.g in fl T , 

(5.1) u n+1 =0 on S T , 
n • Vf +1 =0 on S T , 
u n+1 \ t=0 = m , w" +1 | t=0 = Mi in Q, 
n+1 \t=o = 8o in 0, 

where u n , 6 n , n e N U {0}, are treated as given. 

Moreover, the approximation (u°,9°) is constructed by an extension of the initial 
data in such a way that 

(5.2) M°|t=o=«o, w?|t=o=Mi) 0°\ t= o = 9 in d, 
and 

(5.3) u° = 0, n • V6>° = on S T . 

First we show that the sequence {u n ,9 n } is uniformly bounded. 
Lemma 5.1. (Boundedness of the approximation) Assume that 

D = D(p,p ,q,q ) = \\u \\ W 2 {n} + ||«i|| B 2-2 /P0(fJ) + \\9 \\ B 2-2/ q0(n) 

+ ll&IU PiP0 (n-) + \\g\\L Q , m (n-r) < oo, 
where p,po,q,qo € (l,oo) satisfy the conditions 

3 2 3 2 

- + <1, 

q qo p po 

3 2 3 2 1 

P Po 2q 2(j 

Moreover, assume that 9$> 0_> §, t > § is sufficiently small, and n € NU {0}. T/ien 
£/iere exists a constant A independent of n but depending on D, p, po, q, qo such that 
solutions to problem (5.1) satisfy the estimates 

\\ u t\\wl:l (Q-) + W dn \\w^ (n-) ^ A > 

(5.4) 1 
-e<sup6 n <A 

2 n-r 

for n e NU {0}. 
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Proof. Let u n E W 2 p l o (n T ), 9 n E W^fT). Then Lemma 3.4 ensures the 
existence of solutions to problem (5. 1)1,3,5 an d the estimate 



(5.5) 



r 1 !!^,^) < C (ll W l || W o T ) + ||vr|U P , Po( n T) 

+ ll b IU P , P0 (O-) + l|Wl|l B 2 2/po ), 



where constant c does not depend on r. 
With the use of the formula 



u 



! (x,i) = J v$(x,t')dt' + u (x), 



the first term on the right-hand side of (5.5) is estimated by 



2„.«l 



(5.6) 



\V z u 



J V 2 u?,dt' 



+ r 1 ^||V 2 « || ip (o) 



< r||vVllw^) + t1/p 1I v2w oIIl p( o), 



where in the last line the Holder inequality was used. 

The second term on the right-hand side of (5.5) is estimated with the help of the 
interpolation inequality: 

l|vr|| £piPo(OT) <5||r||^ o( , r) + c (i/5)||r|| L _ (OT) 

(5.7) < %1 } + c(l/8) f 6-dt' + O 

{ JW° T ) 

< ^l^ll^^) +^/S)(T\\O n \\ wlUn , ) +r 1 /«.||eo|U, (n) ) ) 

where c(l/J) ~ <5~ a , a > 0. The interpolation inequality in the first line holds under 
the restriction 



(5.8) 



3 2 3 2 

— 1 < 1, p > q, po>qo, 

q qa v Po 



where the last two restrictions can be relaxed on account of the Holder inequality 
with r < a, ro < er . 

Applying estimates (5.6) and (5.7) in (5.5) leads to 



u 



n+li 



t IIW P ; P0 (O-) < c [ T ll M ril^ :P0 (O-) 

+ (6 + c(l/6)r)\\e n \\ wlUnT) +r 1 /po||V 2 Mo ||i pW 



(5.9) 



c(l/Sy^\\9 Q \\ Lq{n) 



< 4r a (\\u 



af\\ n .n\\ 



-'P.PO 

72,1 



l Ml llB P r P 2 /P0 (n)] 



+ ll«ollw p (Q) + ll M lll B 2 - 2/ ™ 



(Q) 
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where p, p , q, q satisfy (5.8) with p,q,Po,Qo € (l,oo), t is sufficiently small, and 
a > 0. 

By virtue of Lemma 3.3 there exists a solution to problem (5.1)2,4,6, an d the 
following estimate holds 



(5.10) 



\\e 



.71+1 I 



^< V (sup±,sup8")l\\9 n e(u?)\\ Wqo ^ 



+ II |eK)| 2 |U,, TO (f2-) + h\\L q , qo (n-) + Po\\ B i -2 



/q °(ny> 



where ip is an increasing positive function of its arguments and r is sufficiently small. 
By the Holder inequality it follows from (5.10) that 



(5.11) 



\w 



n+l| 



(sup-^sup^W 1 

V fir 9 n OT / 



^2, ,2g 



IVu? 



t lli2,,2, (n T ) 



IVu: 



t \\L 



2q,2q 



(Or) + hh^tl-r) + ||6>0|| S 2 2 /9Q(n) ] 



We proceed now to estimate the norms on the right-hand side of (5.11). First we 
examine 

sup UHlwn) < su P ^||r || B 2-2/, + c(i/<y)||iHk (n) ) 

t<r t<r q ' q a ( ' 



(5.12) 



< sup 



< sup 

t<T 



s\\e n \\ B 2-y qo{n) + c(i/s) 



c(l/S) 




- 8q 











L q (Q). 



t 

5\\6 n || B; -v«o (n) + c(l/S)t^'° ( J (n) df 



1/90 



+ c(l/d)\\e \\ Lg{ n) 
< (S + c(l/Sy^)\\e n \\ w ^ Q{nT) + c (l/J)||6»o|| 2 a 



where c(l/S) <~ S a , a > 0, 1/qo + l/q' = 1. The first inequality in (5.12) expresses 
the interpolation inequality which holds for 

(5.13) 3/q< 2-2/(70. 

The last inequality in (5.12) follows from the imbedding 



(5.14) 



su pII m WIIr 2 - 2 



t<T 



(fi) 



<c( u\ 



ll«(0)|| 



>2-2/ 50 



(n) 



with a constant c independent of r, which holds true for any solution of a parabolic 
equation with vanishing boundary conditions and sufficiently smooth coefficients. 

To estimate 9 n from below by a positive constant we use the assumption 9q > 
9>0. Since by (5.13), 6 n e W%£{W) C C Q (fT) with some a e (0, 1), it follows that 



9 n = 6 n - 9 + 6> > - \9 n -0 \ >£-sup| 



Ic°/ 2 (0,t) t 

Hence for t so small that 

(5-15) \\e n \\ L ^n,c^, T)) r a ' 2 < 49^^^ < U 
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we have 
(5.16) 



> 



To estimate the first two terms in the square bracket on the right-hand 
side of (5.11) we consider the factors ||0 n ||L 2 , >2 , o (^ r ) and ||V<|| £2g2(!o(fr ). 
We have 

iriUwn.) < 5\T\\ W ^) +c(V*)irilw^) 

t 

(5.17) < 5\\e n \\ w , do[nT) + c(l/S) J 6\)dt' + e 



< (5 + c(l/5)r)\\0 n \\ wl , o(nT) +T W°\\0 o \\ Lg(n) 



where the first inequality is the interpolation inequality which holds under the condi- 
tion 

3 2 3 2 

- + --7T-7T- <2, 
q q a 2q 2q 



that is 

(5.18) 

Next, 

(5.19) 



3 2 „ 
- + — < 4. 

q qa 



l|V<IU 2<; , 290 (^) < 5\\un wlUn , } +c(l/J)K|| wnT) 

t 

< S W u tWwlUn n + c(V*) / ul,dt' + ut(O) 



L p,p ( Qt ) 



<(<5 + c(l/*)r)||«?||^ (nT) +c(l/J)r 1 /«.||„ 1 || £p(n) 



where the first inequality holds under the condition 
3 2 3 2 



(5.20) 



p pa 2q 2q 



< 1, p,p ,q,qo e (l,oo). 



Using (5.12), (5.16), (5.17), (5.19) and choosing appropriately S so that to get the 
factors of \\6 n \\ w 2,i ^ QT ^ and ||u™|| w 2,i proportional to r a with a > 0, we infer 
from (5.11) that 



(5.21) 



+ ^H<ll^ o( n T ) + H«olli.(n) + H«illi(o) + IMI W^) 

+ ||6» || B 2 2 /?0(Q) ] 



for q,q ,p,po € (l,oo), restricted by (5.13) and (5.20) ((5.18) is inactive). 
Let us denote 



(5.22) 
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and 

^ D = || M 1| ^(O) + ||ui||^-yPo (n) + ll^|| s 2-2/ 90(n) + ||b|| £ppo( ^) 

+ \\g\\L q , qo (n-y 
Then inequalities (5.9) and (5.21) give 

(5-24) IK^H^,^ < c[r a X n (r) + D], 

and 

(5.25) P n+1 \\wl^) < ^{T a X n {T),D)[{T a X n {r)f+D + D\ 

where p,Po, q 1 <Zo € (l,oo) satisfy (5.8), (5.13) and (5.20), a > 0, r is sufficiently small, 
and c does not depend on r. 

It follows from (5.24) and (5.25) that 

, 26) X n+1 (r) < V (T a X n (r),D)[T a X n (r) + (r a X n (r)) 2 + D + D 2 } 

[ ' ' for n £ NU {0}, 



where 



and p,po, q, qo € (l,oo) satisfy (5.8), (5.13) and (5.20). 

Regarding the fact that <p stands for a generic, positive, increasing function of its 
arguments, (5.26) can be expressed in the following simpler form 

(5.27) X n+1 (r) < ip{T a X n {r),D) for n £ NU {0}. 
For r sufficiently small there exists a positive constant A such that 

(5.28) X (t)<A and ip( T a A,D)<A. 
Then (5.27) implies that 

X n (r) < A for any n £ N. 

This establishes (5.4)i. Moreover, in view of (5.15) and the imbedding W 2 ^ Q (VL T ) C 
C a (Q T ), the bounds (5.4) 2 hold. Thereby the proof is completed. □ 

To show the convergence of the sequence {u n , 9 n } we introduce the differences 

(5.29) U n {t)=u n {t)-u n - 1 {t), d n {t) = e n {t)-6 n - 1 {t), neNU{0}, 
which are solutions to the following problems: 

U r t \ +1 - V • (A l£ (U? +1 )) = V • (A 2 e(U n )) - V • (A 2 mT) in f! T , 

(5.30) U n+1 =0 on S T , 
U n+1 \ t=0 = 0, U? +1 \ t=0 = in fi, 



and 



(5.31) 



c v 9 n ^ +1 - kAd n+1 = -c v d n e n t - 9 n (A 2 a) ■ e{U?) 

- tf n (A 2 a) ■ eiu™- 1 ) + A l£ {U?) • e«) + A^u^ 1 ) ■ e{U?) in ft T , 

n-W n = on S T , 

$ n+1 \t=o = in O, 
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where 

e(u~ 1 ) = 0, e(ut 1 ) = 0, 9^ = 0, 
e(U°) = e(« ), e(t7?) = e(«i), ^° = 6 . 
Let us introduce the quantity 

(5-32) Y"(t) = \\Un wrpo ^) + W^Ww'faar), 

where p,po,q,qa e (l,oo). 

Lemma 5.2. (Convergence of the approximation) Let the assumptions of Lemma 
5.1 hold. Moreover, let the numbers p,po,q,q~o € (l,oo) satisfy the following condi- 
tions: 



q = 2(?, qo = 2q Q 

3 2 3 2 

Q 9o P Po 

3 2 3 2 

P Po 2q 2q 

3 2 3 2 

p po 2q 2q 



< 1, 



< 1, 

< 1. 



TTien there exists a positive constant d which depends on A, D, 0, p,p ,q, 
qo,p,po such that 

(5.33) y™ +1 (r) < dT a Y n ( T ), 
where n <G N U {0}, a > 0, and 

y° = ||n ||^ o(sr) + ||^||^ (OT) . 

Proof. By Lemma 3.4 the solutions to problem (5.30) satisfy the estimate 

(5.34) W^Ww^w < Cl (||V 2 l/l W° T ) + ||W"|| wa .)), 
where c\ does not depend on r. 

The first term on the right-hand side of (5.34) is estimated by 



|V 2 [T' 



V 2 U?,dt' 



<cr\\Un wkk ^y 



and the second one by 

<(* + c(l/<5)r)||<? n ||^ o(nT)) 

where c(l/S) = S~ a , a > 0, and the interpolation inequality in the first line holds 
under the restriction 



(5.35) 



3 2 3 2 

q % P Po 



< L p,Po,q,qo e (i,oo). 
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Hence, 



(5.36) "" p - Po 
< cr a Y n (r). 

For the solutions to problem (5.31) we have 



IK +1 H^rnM<W 



ip sup a , sup ■ 



|tf"0"IU,-, fo (n-) 



(5-37) + ||0»|VI7?| II W^IKI^" 1 ! || wtlT) 

+ || |VC/ri |V<| || W no + II |VE7?| iV^r 1 ! || wir) 

By virtue of Lemma 5.1 the arguments of the function ip are estimated by A and 8/2. 
The successive terms under the square bracket on the right-hand side of (5.37) arc 
estimated as follows. 
The first term by 

II^IU 2f , 2?0 ^)||^IU 2f , 2?0 (n T ) < p(A)r a \\^\\ w ^ y 

This estimate follows on account of (5.4)i and similar arguments as in (5.7) under the 
restrictions 

q = 2q, q = 2q Q , and 

(5.38) 3 2 3 2 3 2 

- + — - xr-^<2, so — + — < 2. 
q qo 2q 2q 2q 2q 

The second term is estimated by 

ll e "IU2,,2 ?0 (n-)l|V(7"|| L2 _ >2 _ o(OT) < ||6»™|| L(!i5o(OT) ■ T a \\Ut\\ w 2,i_ g(nT) 

< V (A)T«\\Un wik{nni 

where we assumed that 

3 2 3 2 

(5.39) - + — < 1. 

p po 2q 2q 

The third term is estimated by 

where we assumed (5.38) and applied (5.4) under the condition 

(5.40) 3 + A_A_JL<i. 

p po 2q 2q 
In view of (5.39) and (5.40) the fourth term is estimated by 

A similar estimate holds for the last fifth term. 
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Using the above estimates in (5.37) yields 
(5.41) \\^ +1 \\w^) < v{A,DA)T a Y n {r). 

^From (5.36) and (5.41) we conclude (5.33). The proof is complete. □ 
Lemmas 5.1 and 5.2 imply the existence of a local solution. 
Lemma 5.3. (Local existence) 

Assume that u G W£(fi), m £ B 2 p - 2 J Po (fl), 6 E B 2 ; q 2/qo (n), 

b e L P:P0 (Q T ), g e L qm {VT), 9 > 6 > 6* > 0, where 9* is given by (4-2), 

p,pa,q,qa G (l,oo) satisfy 

3 2 3 2 3 2 

- + < 1, - + — < 2, 

/e 49] q 9o P Pa P Pa 

{ ' 3 2 3 2 3 2 

- + — - — -—< 1, - + -<2. 
V Pa 2g 2q q q 

Then for r sufficiently small there exists a solution to problem (l.l)-(H) such that 



and 
(5.43) 
where 
(5.44) 



«teW^(n 9€^(n ueC([0,r];^((l)), 



IKHwt >P0 (o-) + H llw, 2 ,'i(n-) + ll«llc([o,r] ; wJ(n)) 
< <p(D,6*,p,p ,q,q ), 



D =D{p,p Q ,q,q Q ,T) = \\u \\ W 2 (n) + \\m || B 2-2/ P0(n) 
+ ll 6 'o|| i j2-2/, 0(! - 2) + ||6|| ipiPo (nr) + ||5||i,, i90 (nT). 



Proof. To satisfy the assumptions of Lemmas 5.1 and 5.2 we choose q = 2q, 
<7o = 2</o and | + > | + Moreover, the assumptions of Lemma 5.1 imply that 
- + — <4, - + — <2, and Lemma 5.2 gives - + — <2, I + J-<1. Hence, we choose 

q go 1 p pa ' ° q qo ' P Po ' 

| + ^ < 2, | + ^ < 1, | + ^ < 2. Then the sequence «,6»"), n <= N, of solutions 
to the approximate problem (5.1) is for sufficiently small r uniformly bounded in the 
space (fT) x (JT ) and convergent in W*^ x ^ (0 T ). 

Hence, by virtue of the well known result (see the proof of Lemma 2.1 from [15, Ch. 
2]) we conclude that the limit 

( Ut = hm ui\e= hm e n ) e w 2 p :l o (n T ) x w, 2 ^(fi T ) 

and satisfies problem (1.1)— (1.4). □ 

Remark 5.4. Under the assumptions 

ueC([0,T];W 2 p (n)), u t e C([0,T] ;J B 2 : p 2 ^(f7)), 

^C([0,r];4 2 /'»(!])), 6ei„(fix(M + T)), 

9 €I M0 (!1x(MIt)), 6{t)>6* >0, te[0,T], 
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with p,Po, q, <7o satisfying (5.42), the assertion of Lemma 5.3 holds true for the interval 
[t,t + t], where t is a sufficiently small number depending on the data. 
Lemma 5.5. (Global existence) Assume that 

u eW 2 r (n), Ul eB 2 r - 2 J r °(Q), 9 eCf°(fi), 
g e L a ^ (n T ) n L oo , r0 (0 T ), b e L r , ro (n T ), 

with r > max{12,p} 7 ro > max{12,po}; °~ > max{6,<7}, ctq > max{6,<7o}, where 
(p,Po) and (q,qo) satisfy (5.42). Then there exists a solution to problem (1.1)-(1.4) 
such that 

Ut ew^. o (n T ), eew 2 'X(n T ), ueC([o,T];W?(n)), 

satisfying the estimate 

^ ll M *llw^ (or) + ll M llc([o/r] ; w^(f2)) + \\ e \\wS : l (nT) 

< tp(D(r, r ,o-,o- ,T),6*,r, r , cr, a , T), 

where D(r, r , cr, CTo, T) is defined by (5.44) an d f is a generic function. 

Proof. To extend the local solution from Lemma 5.3 step by step we choose 
parameters p,po,q,qo in such a way to satisfy simultaneously the restriction (4.115) 
required by a priori estimate (4.116) and the assumptions of Lemma 5.3. Then, in 
accord with (4.115) we assume that 

ll^ollw^ (fi)n^(n) < c \\ u o\\w 2 r (n) 



for r > max t ^ , p , 



(5.46) 

H M lHB^ /™ (o)nB ii/6 2(o) < c||tti|| B;; -v^ (n) 

for r > max{12,p}, r > max{12,p }, 

for a > max{6, q}, ctq > max{6,<7o}- 

For r and a specified in (5.46), the conditions (5.42) are satisfied, because - + — < ^ 
and - + —<§. " 

cr (To — 6 

Let t € [0, T]. By (4.116) in Corollary 4.25 and by the direct trace theorem we 
have 



(5.47) 
Moreover, 



< tp(D,6*,r,r ,o-,o- ,T). 



t 

\\u{t)\\ W l {Qj) < J \\Uf{t')\\ W 2 m dt' + \\uo\\ W 2 {n) 


7 1 1/ 

^ 5 - 48 ) < T i-iAo^| || Ut ,(t')||^ (n) dt') + ||«o||^) 



< T 1 1/ro \\u t \\ w 2,i o(nT) + \\u \\ w i(n) 

< ^(Do,^*,?-,^,^^,? 1 ) 
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for any t € [0, T] . Let N e N be given so large that -^q-j- < r, where t is determined by 
Remark 5.4 and Lemma 5.3. Using (5.47), (5.48) for t = kr and taking into account 
that 

SUP ||6|U (fix(fcT,(fc+l)r)) < \\b\\L riro (nT), 
0<k<N 

SUP ||fl l ||l, . i0 . (nx(fcT,(fc+l)T))nI,oo,ro("x(fer,(fe+l)T)) 

0<k<N 



< 



IMU <r , ao (fi')nz, 00 , ro (n T ) ; 



it follows from Remark 5.4 that there exists a solution (u, 9) on the interval [kr, (k + 
1) T L < fc < A^, where r does not depend on k. 

This implies the existence of the solution of problem (1.1)— (1.4) on the whole interval 
[0,T], which ends the proof. □ 

Choosing r = ro = 12, a = oo = 6 in Lemma 5.5 and recalling Lemma 4.1 we 
complete the proof of Theorem A. 

6. Uniqueness (Proof of Theorem B). Assume that we have two solutions 
(tii,#i), (u 2 ,9 2 ) to problem (1.1)— (1.4). Introducing the differences 

(6.1) U = U!-u 2 , t? = 0i-02, E = e 1 -e 2 , 
we see that they satisfy the problem 

(6.2) U tt -V ■[A 1 E t + A 2 (E-#a)] = 0, 



(6.3) 



c v (6 2 # t + iWit) - kAd = -[0(A 2 a) ■ e lt + 6 2 {A 2 a) ■ E t ] 
+ A\E t ■ en + Ai£ 2t ■ E t in T , 



(6.4) U\ s t=0, n-Wd\ S T=0 on S T , 

(6.5) U\t=o = 0, C/ t |t=o = 0, 0\t=o = O in 0. 
Multiplying (6.2) by U tl integrating over and integrating by parts we obtain 



ld_ ( 
2dt\ 



J U 2 t dx + J (A 2 E) ■ Edx^j + J A\E t ■ E t dx = j {A 2 da) ■ E t dx. 



n n 
In view of (1.10) this implies 



\j t (y + y ( A 2-E) • Edx^j + a u y |-E t | 2 cfe < y (A 2 ??o;) • B t da: 
an si si 

which, after applying the Holder and the Young inequalities to the right hand side, 
yields 

(6.6) ~ y (U 2 + (A 2 E) ■ E)dx + ^a u j \E t \ 2 dx < a j d 2 dx. 

n n q 



(6.7) 

+ c 2 / (^\e lt \+0 2 \E t \\m)dx 
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Further, multiplying (6.3) by integrating over fi, integrating by parts and using 
the boundary conditions implies 

\ j c v (e 2 ^$ 2 + 8 lt V 2 )dx + k j\sJd\ 2 dx 

n n 

<c 2 J{d 2 \e lt \+0 2 \E t \\$\)dx. 
n 

Continuing, we have 

Ijt J c - e ^ 2dx + k J \W\ 2 dx<c 3 j{\9 lt \ + \8 2t \)d 2 dx 
inn 

c 2 J(ti 2 \e lt \+8 2 \E t \\d\)da 
n 

^From (6.7) we get 

\jt I c ^ 2dx + k Mm(n)<^\\nl 6i n) 
n 

(6.8) + c(l/ei)(||M£ 3 (n) + IIM£: 

+ e2||^||L(o)+c(l/£2)||eit||l3(o)||^HL(o) 
+ £3\\E t \\ 2 L2{n) + c(l/ £ 3)||^||L(n)ll^llL(o) + m\l 2{ ny 

Adding (6.6) and (6.8), assuming that s\, e 2l £3 are sufficiently small we obtain 

j t J(U 2 + (A 2 E) ■ E + c v 6 2 $ 2 )dx + ^a u J \E t \ 2 dx 

(6 9) n n 

+ W\\ 2 HH n) < c(ci + ||M£, ( n) + ||fl«||i 3(n) + ||ei t ||i 3(n) 

+ IMl^n) +k)\\nl 2{ ny 

By virtue of (4.2) it holds 

6» 2 > 0, > 0. 

Thus, introducing 

X{t) = j(U 2 t + (A 2 E) ■ E + c v 9 2 d 2 )dx, 
n 

A(t) = Cl + k + ||Mis(n) + IIM£ 3 (n) + lkit||| 3 (Q) + IMlLtn). 
we conclude from (6.9) the inequality 

(6.10) j t X<AX for te (0,T). 
Hence, 

t 

(6.11) X(t) < X(0) cxp ^ J A{t')dt^j . 
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Since X(0) = and, by the assumption (1.17), J A(t')dt' < oo, it follows that 
X(t) = for t <G (0,T). This proves the uniquenes of the solution. 
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